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ABSTRACT 

The radio counterparts to the IRAS Redshift Survey galaxies are identified in the NRAO 
VLA Sky Survey (NVSS) catalog. Our new catalog of the infrared flux-limited (S'go/jm > 2 Jy) 
complete sample of 1809 galaxies lists accurate radio positions, redshifts, and 1.4 GHz radio and 
IRAS fluxes. This sample is six times larger in size and five times deeper in redshift coverage 
(to z K, 0.15) compared with those used in earlier studies of the radio and far-infrared (FIR) 
properties of galaxies in the local volume. The well known radio-FIR correlation is obeyed 
by the overwhelming majority (> 98%) of the infrared-selected galaxies, and the radio AGNs 
identified by their excess radio emission constitute only about 1% of the sample, independent of 
ttie IR luminosity. These FIR-selected galaxies can account for the entire population of late-type 
field galaxies in the local volume, and their radio continuum may be used directly to infer the 
extinction-free star formation rate in most cases. Both the 1.4 GHz radio and 60 /im infrared 
luminosity functions are reasonably well described by linear sums of two Schechter functions, 
one representing normal, late-type field galaxies and the second representing starbursts and 
other luminous infrared galaxies. The integrated FIR luminosity density for the local volume 
is (4.8 ± 0.5) X IO^Lq Mpc^'^, less than 10% of which is contributed by the luminous infrared 
galaxies with LpiB. ^ IO^^Lq. The inferred extinction-free star formation density for the local 
volume is 0.015 ± 0.005 Mq yr^^ Mpc'^. 

Subject headings: galaxies: luminosity function - galaxies: starburst - galaxies: active - infrared: 
galaxies - radio continuum: galaxies — surveys 



1. Introduction 

The Infrared Astronomical Satellite (IRAS; 
Neugebauer et al. 1984) was the first telescope 
with sufficient sensitivity to detect a large number 
of extragalactic sources at mid- and far-infrared 
wavelengths. The IRAS survey covering 96% of 
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the sky detected a sample of ~ 20, 000 galaxies 
complete to 0.5 Jy'^ in the 60 /xm band, the ma- 
jority of which had not been previously cataloged. 
While most of these are late-type galaxies with 
modest infrared luminosity {Liu/Lb = 0.1 — 5), 
some extremely bright galaxies emitting most of 
their light in the infrared wavelengths were also 
detected {Ljr > 10^^ Lq, Liu/Lb > 50; de Jong 
et al. 1984, Soifer et al. 1984, 1987 - see review by 
Sanders & Mirabel 1996). Their observed FIR lu- 
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minosity can reasonably be explained by ongoing 
star formation and starburst activity, but some of 
the most luminous infrared galaxies may host an 
active galactic nucleus (AGN). 

There are several important reasons for investi- 
gating the radio properties of FIR-emitting galax- 
ies. First of all, the observed radio emission may 
be used as a direct probe of the very recent star 
forming activity in normal and starburst galax- 
ies (see a review by Condon 1992). Nearly all of 
the radio emission at wavelengths longer than a 
few cm from such galaxies is synchrotron radia- 
tion from rclativistic electrons and free-free emis- 
sion from H II regions. Only massive stars (M ^ 
8 Mq) produce the Type II and Type lb super- 
novae whose remnants (SNRs) are thought to ac- 
celerate most of the relativistic electrons, and the 
same massive stars also dominate the ionization 
in H II regions. This results in an extremely good 
correlation between radio and far-infrared emis- 
sion among galaxies (Harwit & Pacini 1975, Con- 
don et al. 1991a, b). While the use of radio contin- 
uum as a star-formation tracer has been suggested 
before, the validity and details of the method has 
not yet been quantified fully. By examining the 
radio properties of a large sample of FIR-selected 
galaxies, we plan to address the utility of radio 
continuum emission as a probe of star-forming ac- 
tivity in detail, including its limitations and the 
magnitude of systematic errors. 

Another important reason for examining their 
radio properties is to investigate the frequency 
and energetic importance of AGNs among FIR- 
bright galaxies. If a luminous AGN is present, 
then the associated radio continuum, in excess of 
the level expected from the star formation, may be 
detected. The presence of a moderately compact 
nuclear radio continuum source is not a sufficient 
proof that an AGN is the source of the observed 
luminosity. On the other hand, absence of signifi- 
cant radio emission associated with AGN activity 
may indicate that AGNs are not energetically im- 
portant, at least among the FIR-selected galaxies 
(e.g. Crawford ct al. 1996). 

This radio study of the FIR-selected galaxies is 
particularly timely and relevant in light of the re- 
cent exciting discovery of highly obscured and lu- 
minous submm sources at cosmological distances 
(Small, Ivison, & Blain 1997; Hughes et al. 1998; 
Barger et al. 1999; Eales et al. 1999). Though 



optically unremarkable, these dusty galaxies may 
dominate the star formation and metal production 
at ^; > 2. The majority of these submm galax- 
ies may be too faint optically to yield their red- 
shifts even with the largest telescopes, and deep 
radio imaging observations such as with the VLA 
may be the only way of identifying and inferring 
their redshift distribution and star formation rate 
for now (CarilU & Yun 1999, Small et al. 2000). 
Therefore, understanding the radio properties of 
the FIR/submm-selected galaxies in the local uni- 
verse may offer a valuable insight into interpreting 
the nature of these highly obscured galaxies in the 
early universe. 

In the study of an optically selected sample of 
299 normal and irregular galaxies, Condon, An- 
derson, & Helou (1991a) found a non-linear rela- 
tion between radio, FIR, and B-band luminosity, 
independent of Hubble type. Studies of an IRAS- 
selected galaxy sample by Lawrence et al. (1986; 
303 galaxies in the North Galactic Pole), Unger 
et al. (1989; 156 galaxies in the South Galactic 
Pole), and Condon et al. (1991a; 313 galaxies in 
the IRAS Bright Galaxy Sample [BGS]) extended 
the observed tight correlation between radio and 
FIR luminosity to the most luminous galaxies in 
the local universe. The FIR emission of galaxies is 
sensitive to environmental effects, and density in- 
homogeneity such as clustering can complicate the 
derivation of a luminosity function (e.g. Oliver et 
al. 1996). Therefore a study involving the largest 
possible volume, bigger than large-scale structures 
such as cosmic walls and voids, is highly desirable. 
The aim of this paper is to conduct a new exten- 
sive analysis of radio and FIR emission in star- 
forming galaxies using the largest sample avail- 
able. Our IRAS-selected sample with full redshift 
data is six times larger (1809 galaxies) in size and 
five times deeper in distance than all similar stud- 
ies and covers over 60% of the sky. Therefore our 
results should be more immune to such spatial in- 
homogeneity and statistical fluctuations. Both the 
normal late- type galaxies and interaction/mcrgcr- 
induced high luminosity galaxies are included in 
our sample, and possible differences between their 
luminosity functions are examined. 

An accurate derivation of the local radio and 
FIR luminosity functions also allows the determi- 
nation of the local radio and FIR luminosity den- 
sities and thus an extinction-free estimate of the 
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local star formation rate (SFR). For this reason, 
all cluster galaxies are also explicitly included in 
our analysis, unlike some of the previous studies. 
Most existing estimates of the local SFR are based 
on either UV continuum or optical emission line 
studies (Gallego et al. 1995; Lilly et al. 1996; 
Treyer et al. 1998; Tresse & Maddox 1998) that 
are subject to highly uncertain extinction correc- 
tions (e.g. Meurer et al. 1997, Steidel et al. 1999). 
Optically thin FIR and radio continuum emission 
provides an inherently extinction-free estimate of 
star formation rate and a good calibration for the 
UV and optical techniques. 

The organization of the paper is as follows. The 
sample selection and the correlation of the radio 
and FIR properties are described in §2. The well 
known radio-FIR correlation is examined in de- 
tail, and "radio-loud," "radio-excess," as well as 
"infrared-excess" galaxies are identified in §3 &: 
§4. The radio and FIR luminosity functions for 
the local volume are derived, and inferences are 
made on the nature of these luminosity functions 
in §5. Finally, the frequency of the AGNs and their 
energetic importance are discussed, and the in- 
frared luminosity density is derived and discussed 
in terms of the star formation rate in the local 
volume in §6. The strengths and limitations of 
using the radio-continuum luminosity density for 
measuring the star formation rate and cosmic star 
formation history are also discussed. 

2. Sample Selection and Data Analysis 

The main catalog used for our sample selec- 
tion is the 1.2 Jy IRAS Redshift Survey catalog'* 
(Strauss et al. 1990, 1992, Fisher et al. 1995). 
This catalog contains 9897 sources selected from 
the IRAS Point Source Catalog (PSC) by the 
following criteria: S'eo^m > 1.2 Jy, (Seo^im)^ > 
(<S'i2,im X 5'25^m), and |6| > 5° (Strauss et al. 
1992). Of the 9897 sources covering 87.6% of the 
sky, approximately 5321 are galaxies and 14 were 
of unknown type as of 1995 February 6. Among 
the non-galaxian sources are HII regions, stars, cir- 
rus, and planetary nebulae. 

^Thc 1.2 Jy IRAS Redshift Survey catalog was obtained 
from the Astronomical Data Center (ADC) operated by 
NASA/Goddard Space Flight Center. 



The matching radio data for the FIR-selected 
sample were obtained from the recently completed 
NRAO VLA Sky Survey (Condon et al. 1998). 
The NRAO VLA Sky Survey was made with D 
and DnC configurations of the Very Large Array at 
a frequency of 1.4 GHz and an angular resolution 
of 45". Observations began in 1993 and extended 
into 1997, covering the celestial sphere above 6 = 
-40°. The final NVSS catalog contains « 1.8 x 10*^ 
sources, including most of the galaxies in the IRAS 
Faint Source Catalog, which has a flux limit of 
Seoiim ~ 0.3 Jy. 

Both the IRAS redshift catalog and the NVSS 
catalog have well-known limitations in their com- 
pleteness (see above). To construct an unbiased 
complete sample for statistical analysis, we se- 
lected objects at (5 > —40° and located away 
from the Galactic plane (|6| > 10°). The red- 
shift data are complete only for the subsample 
with S'eo/am > 2 Jy (Strauss et al. 1992), and 
we adopted this flux cutoff for completeness. The 
sensitivity of the NVSS survey is about 0.5 mJy 
beam~^ {Icr), sufficient to detect all sample galax- 
ies with S/N > 20. 

We first constructed a catalog of FIR-selected 

galaxies with 1.4 GHz radio counterparts by cross- 
correlating the 9897 sources in the 1.2 Jy IRAS 
Redshift Survey with the radio sources in the 
NVSS database using a program based on a soft- 
ware written by W. Cotton. We sought radio 
counterparts in the NVSS catalog within 30" of 
the IRAS PSC positions, and unique counterparts 
were found in nearly all eases. This search radius 
is twice the size of the typical la error ellipse in 
the IRAS PSC, and the probability of finding an 
unrelated NVSS source within each search radius 
is less than 6 x 10^'' (see Langston et al. 1990). 
Two or more radio counterparts were found in 19 
cases - about three times larger than expected 
from chance superposition. Examination of these 
multiple matches revealed that they are mostly 
multiple NVSS components produced by a single 
galaxy with a large angular extent (e.g., M82). 
Because the NVSS is a snapshot survey, some 
of the extended disk emission in galaxies with a 
large angular extent was missed. Therefore the 
1.4 GHz fluxes for the 134 large late-type galaxies 
and multiple matches were replaced with the mea- 
surements made in the multiple-snapshot study of 
bright spiral galaxies and the IRAS BGS galaxies 
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by Condon (1987) and Condon ct al. (1990). Only 
25 cases required corrections larger than 30%. The 
IRAS fluxes are from the PSC, but those for the 
galaxies with optical diameter larger than 8' are 
replaced by the imaging measurements by Rice et 
al. (1988). 

The 1.4 GHz and 60 fim luminosities were cal- 
culated using the following relations: 

log Li iGHzi^ Hz"^) = 20.08+2/os'£)+/os'S'i 4ghz 

(1) 

log (Lq) = 6.014 + 2logD + logSeo^m (2) 

where D is the distance in Mpc^, and S\aghz and 
SeOfim are flux densities in units of Jy. This 60 /im 
luminosity (ieo/nm) is the luminosity contribution 
from the IRAS 60 /xm band to the FIR luminosity 
Lfir^, i.e. 



{Lq) = (1 + 



2.58S'60/it 



(3) 



(sec Hclou et al. 1988) and is related to an- 
other commonly cited quantity "i/L,y(60 /xm)" 
(e.g. Soifer et al. 1987, Saunders et al. 1990) 
as Leopm = 0.69 x vL^{QQ fim). 

Our final catalog of the IRAS 2 Jy Galaxy Sam- 
ple sources meeting our selection criteria contains 
1809 galaxies. The names, radial velocities, opti- 
cal magnitudes, 1.4 GHz and 60 ^m luminosities, 
and FIR/radio flux ratios q (sec Eq. 5) of the 161 
luminous FIR galaxies having Lgo/^ > IO^^'^^Lq 
are given in Table 1. We plan to publish the full 
catalog in an electronic form elsewhere. 

The distribution of 1.4 GHz radio luminosity 
for the IRAS 2 Jy sample is shown as a function 
of galaxy redshifts in Figure 2. Low-luminosity 
galaxies drop out quickly with increasing redshift 
due to the IRAS 60 /im flux cutoff at 2 Jy, but the 
infrared luminous galaxies are present to z = 0.16. 
Therefore, our IRAS 2 Jy sample is not only 6 
times larger in the sample size compared with 
similar previous studies, it also includes sources 
5 times further away in the redshift domain as 
well. The bona fide radio-loud objects in the sam- 
ple arc PKS 1345-^12, 3C 273, and NGC 1275, 



^Ho = 75 km s ^ Mpc ^ is assumed throughout this paper. 

^The IRAS 100 fiin flux is not available for a small subset 
of the IRAS redshift catalog sources, but all such sources 
axe excluded from our complete sample by other selection 
criteria. 



all of which stand out clearly in Figure 2 with 
LiAGHz > 1-5 X 10^5 W Hz-^ Determining the 
fraction of infrared-luminous objects that are also 
radio loud (due to the presence of an AGN) is 
one of the key questions we address later (see §4) . 
One immediate result from this study is the de- 
termination of accurate positions for the IRAS 
Redshift Survey sources. The positions of IRAS 
sources such as found in the IRAS Point Source 
Catalog arc only accurate to about 15" (Ic). In 
comparison, the source coordinates given by the 
NVSS catalog have ~ 1" rms uncertainties, suffi- 
cient for a secure identification of the sources even 
in a crowded or confused field. For example, the 
IRAS position error ellipse for the interacting pair 
IRAS 12173-3514 falls between the two widely 
separated optical galaxies (see Figure 3) while the 
NVSS catalog identifies the northern galaxy as the 
dominant source. 

The histogram of position offsets between the 
IRAS PSC and NVSS (Figure 4) shows that the 
position offset is less than 10" in over 2/3 of all 
cases, and this is consistent with the expectation 
from the IRAS position error. 

3. Radio-FIR Correlation 

The correlation between the radio and FIR lu- 
minosity in galaxies is one of the tightest and 
best-studied in astrophysics. It holds over five or- 
ders of magnitude in luminosity (Price & Duric 
1992), and its nearly linear nature is interpreted 
as a direct relationship between star formation and 
cosmic-ray production (Harwit & Pacini 1975; 
Rickard & Harvey 1984; de Jong et al. 1985; 
Helou et al. 1985; Wunderlich & Klein 1988). 

Evidence for non-linearity in this relation was 
discussed by Fitt, Alexander, & Cox (1988) and 
Cox et al. (1988) for magnitude-limited samples of 
spiral galaxies. Price & Duric (1992) later showed 
that radio emission can be separated into free-free 
and synchrotron emission and suggested that the 
latter causes a greater than unity slope at lower 
frequencies. 

3.1. Observed Radio-FIR Correlation 

The plot of 1.4 GHz radio luminosity versus 
60 /Ltm FIR luminosity for our IRAS 2 Jy sample 
galaxies (Figure 5a) shows a nearly linear correla- 
tion spanning over 5 decades in luminosity. The 
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three noteworthy features are: (1) a unity slope; 
(2) steepening of the relation at low luminosity 
{Leoiirn < IO^Lq); and (3) a higher dispersion for 
huninosities above Leofim ^ IQ-'^'^'^Lq. 

A formal fit to the observed radio-FIR luminos- 
ity correlation yields 

log{LiAGHz) = [0.99±0.01]Zo5(L60/.m/io)+[12.07±0.08]. 

(4) 

The overall trend is indistinguishable from a pure 
linear relation since the overwhelming major- 
ity of FIR-selected galaxies have L^o^m between 
-[^09-11.5^^ The observed strong correlation is 
accentuated by the distance effect, but the plot of 
the observed radio and 60 jim fluxes (Figure 5b) 
also shows a linear trend spanning nearly three 
orders of magnitude in flux. In both plots, the 
scatter is about 0.26 in dex. 

There is a systematic tendency of galaxies with 
Leoum < 10^ to appear below the best-flt line 
in Figure 5a, and the diminished radio emission 
of these galaxies may contribute to the previ- 
ously reported "non-linear" trend in the radio- 
FIR relation among optically selected galaxy sam- 
ples (e.g. Condon ct al. 1991a). Such a de- 
viation from the linear relation can occur if the 
FIR or radio luminosity is not directly propor- 
tional to the star forming activity. Dust heating by 
low mass stars ( "cirrus" emission) may contribute 
only to the observed FIR emission (Helou 1986; 
Lonsdale-Persson & Helou 1987; Fitt, Alexan- 
der, & Cox 1988; Cox et al. 1988; Devereux 
& Eales 1989). Alternatively, the radio luminos- 
ity of less luminous (and generally less massive) 
galaxies may be low if cosmic-ray loss by diffusion 
is important (Klein et al. 1984; Chi & Wolfendale 
1990). Because this sample is FIR selected, there 
is a potential bias towards sources with higher 
FIR/radio flux density ratios (see Condon & Brod- 
erick 1986). On the other hand, low FIR luminos- 
ity sources with relatively high FIR/radio ratios 
may be missed systematically, possibly offsetting 
the first potential bias. 

3.2. Deviation from the Linear Relation 

The presence of any non-linearity or any in- 
crease in the dispersion of the radio-FIR corre- 
lation may be seen more easily by examining the 



"g" parameter (Condon et al. 1991a), 



q = log[ 



FIR 



3.75 X 1012W m" 



log[ 



Si. 



4GHz 



W m-^Hz^i 



(5) 

where Siaghz is the observed 1.4 GHz flux density 
in units of W m~^ Hz"^ and 

FIR = 1.26X 10-i^(2.5856o^„ + 5ioo^„) W m'^, 

(6) 

where /Seo^jm and S'loojum are IRAS 60 /um and 
100 /um band flux density in Jy (see Helou ct al. 
1988). Therefore g is a measure of the logarith- 
mic FIR/radio flux-density ratio. For most galax- 
ies in the IRAS Bright Galaxy Sample q « 2.35, 
but some galaxies have smaller q values due to 
an additional contribution from compact radio 
cores and radio jets/lobes (Sanders & Mirabel 
1996). Optically selected starburst galaxies (inter- 
acting/Markarian galaxies) have about the same 
q values as normal galaxies, and the radio-FIR 
flux ratio appears to be independent of starburst 
strength (e.g. Lisenfeld, Volk, & Xu 1996). 

The q values for the IRAS 2 Jy sample were 
computed using Eq. 5 and are shown in Figure 6. 
The mean q value for the entire sample is 2.34 ± 
0.01, which is in good agreement with q values 
typically found in other star-forming galaxies. In 
fact, the 98% of galaxies in our IRAS 2 Jy sample 
are found between the two dotted lines which mark 
the five times excess and deficit of radio emission 
with respect to mean q = 2.34 (solid horizontal 
line). 

The effects of measurement errors and source 
confusion can be estimated from the dispersion in 
the q value in Figure 6. Before the radio fluxes for 
galaxies with large angular extent were corrected 
using the measurements by Condon (1987), several 
objects with unusually large q values (g > 3) were 
found - resulting from the missing extended flux. 
One object, NGC 5195, also turned out to include 
incorrect IRAS fluxes in the IRAS redshift cat- 
alog, confused with its brighter companion M51. 
Using the correct IRAS fluxes, its g = 2.68 became 
much closer to the sample mean value. Similarly, 
IRAS 02483-1-4302 is a well known quasar-galaxy 
pair and another case of source confusion. Its 1.4 
GHz flux was corrected using the measurement 
by Crawford et al. (1996). These occurrences 
are consistent with the expected source confusion 
statistics discussed earlier (§2.3). The direct com- 
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parison of Figure 6 with a similar plot in Condon 
et al. (1991a) suggests that the measurement er- 
rors and confusion effects are nevertheless small 
and comparable to those of the bright spiral sam- 
ple and the IRAS BGS sample analyzed by Con- 
don et al. 

There are good reasons to believe that galaxies 
with much larger or smaller q values are inher- 
ently different. Firstly, there are several objects 
whose smaller q values are the result of signifi- 
cant excess radio emission associated with AGNs 
(sec §4 for a further discussion). There are also 
a few infrared-excess objects which may be highly 
obscured compact starbursts or dust enshrouded 
AGNs (see §6.1 & §6.2), and some of the observed 
dispersion may arise at least in part from the vari- 
ation in excitation or dispersion in dust temper- 
ature. The clustering of low-luminosity galaxies 
{L^otmi ~ 10^ ^0) above the mean q = 2.34 in 
Figure 6 is another clue that the scatter in the 
q distribution is more than simply statistical in 
nature. This stems from the previously noted de- 
viation from the linear radio-FIR relation (§3.1), 
and it is also the source of a small gradient in 
q noted by Condon et al., due to a systematic 
decrease in radio emission among low luminosity 
galaxies. Lastly, the dispersion in the q values 
among the objects with ieopm > IO^^Lq is sig- 
nificantly larger, ag = 0.33. This larger scatter 
also manifests itself in Figure 6 as the disappear- 
ance of the dense core in the distribution near the 
mean. Helou et al. (1985) noted a similar increase 
in dispersion at high luminosity in their analysis 
of only 38 spiral galaxies and suggested that the 
radio-FIR correlation may break down for high lu- 
minosity sources. The new data suggest that the 
correlation still holds but the dispersion is indeed 
increased (see below for possible explanations). 

4. Radio-Loud and Radio-Excess GalsLxies 

Two different functional definitions of "radio- 
loud" objects are found in the literature: one 
based on absolute radio power and the second 
based on a relative flux ratio (e.g. Lradiol Lopt)- 
The former are generally a subset of the latter. 
Here we call "radio-loud" only those galaxies host- 
ing radio sources with Li^ghz > 10^^ W Hz~^. 
In contrast, we define "radio-excess" objects as 
galaxies whose radio luminosity is five times the 



value predicted by the radio-FIR correlation or 
larger (below the lower dotted line in Fig. 6, i.e. 
q < 1.64). 

4.1. Radio-Loud Objects 

It is clear from Figure 2 that there are only 
three radio-loud objects in our sample of 1809 
galaxies, with several others approaching the 
radio-loud limit. Regardless of the exact statistics, 
the bona fide radio-loud galaxies are extremely rare 
among the infrared selected sample. We discuss 
below the three radio-loud objects in our original 
sample in detail. It appears that each object rep- 
resents a special case, and their radio and FIR 
emission have physically distinct origins. 
PKS 1226-1-023 (3C 273) is a well-known QSO 
at z 0.158 (Schmidt 1963). It is a warm IRAS 
source (S'eo/im/S'ioo^im = 0.76, Tdust ^ 44 K) 
and also the most luminous (log Li,4ghz = 27.42 
W Hz""'^) radio source in our sample. Its ellip- 
tical host is somewhat brighter than the bright- 
est galaxy in a rich cluster (My = —22.1, Bah- 
call et al. 1997) and is probably located inside a 
poor cluster (Stockton 1980). The large-scale ra- 
dio structure from the VLA and MERLIN shows 
a compact, flat-spcctrum core and a single jet ex- 
tending about 23" from the core at a position 
angle of 222° (Conway et al. 1993). Unlike the 
other two radio-loud objects discussed below, its 
infrared emission appears to be mostly the contin- 
uation of the non-thermal spectrum from optical 
to radio wavelengths, associated directly with the 
AGN. 

NGC 1275 (3C 84) is the central giant ellipti- 
cal galaxy of the Perseus Cluster with a Seyfert 2 
nucleus and may be accreting gas from the X-ray 
emitting intercluster medium via a cooling flow 
(Kent & Sargent 1979; Fabian et al. 1981). The 
bulk of its 1.4 GHz radio continuum emission (log 
LiAGHz = 25.13 W Hz~^) comes from the core 
component less than 1" in size, but it also has a 
10' size extended radio structure that may be in- 
terpreted as an asymmetrical Fanaroff-Riley type 
I source whose jet axis lies close to the line-of-sight 
(Miley & Perola 1975; Pedlar et al. 1990). The 
FIR emission is attributed to an extended (^ 30") 
distribution of dust, possibly from a kpc scale 
circumnuclear gas/dust structure accreted from a 
cluster interloper (Lester et al. 1995; Inoue et al. 
1996). Unlike the compact nuclear radio source. 
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the 100 /im emission in NGC 1275 shows httle ev- 
idence for variabiUty (Lester et al. 1995). There- 
fore, the FIR emission in NGC 1275 may be un- 
related to the radio emission, which is dominated 
by the non-thermal emission from the AGN. 
PKS 1345-1-12 (4C 12.50) is the second 
most luminous radio source in our sample (log 
LiAGHz = 26.18). The optical galaxy has a dou- 
ble nucleus embedded in a distorted halo ~ 85 kpc 
in diameter, indicating a merger system approach- 
ing the appearance of a cD galaxy (Mirabel et al. 
1989; Shaw, Tzioumis, & Pedlar 1992). This 
galaxy has a core-halo radio structure, and more 
than 95% of the flux at 2 cm wavelength is com- 
ing from a region smaller than 0".05 (106 pc; van 
Breugel et al. 1984). The western nucleus shows 
a Seyfert 2 spectrum (Gilmore & Shaw 1986), 
and Veilleux et al. (1997) reported broad infrared 
recombination lines indicative of a hidden QSO. 
The eastern nucleus is associated with a compact 
steep-spectrum radio source with a double-lobe 
structiue (Shaw et al. 1992). The large infrared 
luminosity (log Lfir{Lq) = 11.9) makes it a 
genuine ultraluminous infrared galaxy with an ex- 
tremely large gas content inferred from the CO 
observations {Mr^ = 6.5 x 10^° M©; Mirabel 
et al. 1989, Evans et al 1999). Therefore it is 
one of the best candidates for studying a possi- 
ble common origin for ultraluminous galaxies and 
powerful radio galaxies. 

4.2. Radio-Excess Objects 

One of the key objectives of this study is a quan- 
titative analysis of the frequency and energetic im- 
portance of AGNs in FIR-selectcd galaxies. As- 
suming that radio excess is an AGN indicator, a 
total of 23 potential AGN hosts with q < 1.64 
are identified in our sample (see Table 2). The 
low frequency (23/1809 = 1.3%) of these radio- 
excess objects indicates that radio AGNs occur 
quite rarely among the infrared selected sample 
of normal and starburst galaxies. Among the 
FIR-luminous galaxies with log Lpm > 11. ^Lq 
listed in Table 1, there are only three radio-excess 
galaxies (IRAS 12173-3541, IRAS 12265-^0219, 
IRAS 13451-1-1232). This fraction, 3/161 = 1.9%, 
is only marginally larger than the radio-excess 
fraction among the entire IRAS selected sample. 
Therefore the radio- excess fraction is generally 
only a few percent and is nearly independent of 



the infrared luminosity. 

How does this result compare with the anal- 
ysis of optically selected galaxy samples? Con- 
don & Broderick (1988) have identified a total of 
176 radio sources with S\aghz > 150 mJy among 
galaxies in Uppsala General Catalogue of Galaxies 
(Nilson 1973; hereafter UGC). Judging by the ra- 
dio morphology, radio-to-infrared flux density ra- 
tio, and infrared spectral index, 64 are identified 
as "starbursts" and 103 (59%) are identified as 
"monsters" . There are four borderline cases while 
no infrared data are available for the remaining 
three. We find that the same distinctions could 
have been made using their radio-to-FIR flux ratio 
alone. Among the "monsters" in the UGC sam- 
ple, only 6 galaxies (UGC 2188, 2669, 3374, 3426, 
8850, and 12608) are bright enough at 60 /xm to 
be included in our sample, and all but UGC 2188 
are identified as radio-excess by our analysis. The 
starburst/Seyfert galaxy UGC 2188 (NGC 1068) 
has q = 1.72 and is also marginally radio-excess 
(about 4 times larger radio luminosity than ex- 
pected from the radio-FIR correlation) . Therefore 
there is little overlap between the "monsters" in 
the UGC sample and our IRAS 2 .ly sample. This 
is mainly due to the relatively large radio flux cut- 
off used in the UGC study, which has introduced a 
strong bias in favor of radio bright AGNs. As the 
radio flux cutoff is lowered, the dominance of the 
star-forming galaxies emerges quickly as shown by 
the preliminary analysis of the UGC galaxies us- 
ing the NVSS database (Cotton & Condon 1998). 
The relative frequencies of "monsters" and star- 
forming galaxies is addressed further in terms of 
their luminosity functions below. 

The frequency of radio-excess objects derived 
above is strictly a lower limit to the overall fre- 
quency of AGNs among the IR selected galax- 
ies. Radio luminosity of AGNs ranges over 10 
orders of magnitudes (e.g. Condon & Broderick 
1988, Ho & Ulvestad 2000), and radio-quiet AGNs 
might not be that rare in our sample. Neverthe- 
less a large majority of luminous AGNs are radio 
sources and may be readily identified. For exam- 
ple, 84% (96/114) of Palomar Bright Quasar Sur- 
vey sources have been detected by Kellermann et 
al. (1989) at 5 GHz using the VLA. Evidence for 
radio emission in excess of the inferred star for- 
mation activity level among high redshift dusty, 
radio-quiet QSOs has also been discussed by Yun 
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et al. (2000). To quantify the robustness of radio- 
excess as an AGN indicator, we have examined the 
radio properties of a sample of 48 Palomar- Green 
(PG) survey QSOs with IRAS detections reported 
by Sanders et al. (1989). Ten of these sources were 
dropped from our analysis, either because the XS- 
CANPI analysis revealed the IRAS detections to 
be the results of confusion or because meaning- 
ful upper limits for the FIR and radio luminosity 
could not be obtained. Among the remaining 38 
PG QSOs, only 23 have Leo^m > 10"i©. The 
radio-excess fraction among this subsample of FIR 
luminous QSOs is 12/23=52%, which greatly ex- 
ceeds the frequency for our IRAS selected galaxy 
sample (see above).'' Among the remaining 11 PG 
QSOs without the excess radio emission, AUoin 
et al. (1992) detected CO emission with inferred 
molecular gas mass larger than lO^^M© in all three 
objects observed and argued that star formation 
may be responsible for the bulk of the FIR lumi- 
nosity among these QSOs. Therefore the overall 
frequency of AGNs may be somewhat larger than 
inferred from the radio-excess alone - perhaps as 
much as twice as large in the case of the IRAS 
selected PG QSOs. Further discussions of the fre- 
quency of energetically dominant AGNs are found 
in the Discussion section (§6), including the exam- 
ination of the mid-infrared properties. 

5. Radio and Far-infrared Luminosity 
Functions 

With 1809 galaxies in our IRAS 60 /xm flux- 
limited sample, we should be able to derive the 
radio and FIR luminosity functions (LFs) for field 
galaxies in the local volume with significantly bet- 
ter statistics than previously possible. Since our 
sample covers a significantly larger solid angle and 
volume than the samples used previously, fluctua- 
tions produced by large-scale structures and clus- 
tering should be smaller. If there is a significant 
presence of radio AGNs contributing to our IRAS 
selected sample, not only they can be identified as 
a separate population in the derived radio and/or 
FIR luminosity functions, we can also determine 
the range of luminosity where their contribution 
is significant. 

The radio-excess fraction for the entire range of LeOMm is 
between 34% and 45%, depending on how the sources with 
only upper limits of q are accounted. 



Several different functional forms have been 
suggested previously for the radio and FIR lu- 
minosity functions of field galaxies, such as log- 
normal (Hummel 1981; Isobe & Feigelson 1992), 
hyperbolic (Condon 1984), two power-laws (Soifer 
et al. 1987), and power law/Gaussian (Sandage 
et al. 1979; Saunders et al. 1990). These func- 
tional forms offer reasonable descriptions of the 
observed radio and FIR LFs with varying degrees 
of relative merit. All previously published stud- 
ies, however, have unanimously rejected the tradi- 
tional form proposed by Schechter (1976) as being 
too narrow to describe the observed radio or FIR 
LFs. This is surprising because Schechter func- 
tion offers an excellent and robust description of 
the optical LFs of field galaxies and galaxy clusters 
(e.g. Shapiro 1971, Schechter 1976, Auriemma et 
al. 1977, Efstathiou, Ellis, & Peterson 1988). Ex- 
tinction effects at optical wavelengths cannot ex- 
plain this discrepancy because Schechter function 
also offers an excellent description of near-IR LFs 
of galaxies (Gardner et al. 1997; Szokoly et al. 
1998; Loveday 2000). Therefore the failure of the 
Schechter form is highly unusual and deserves fur- 
ther investigation. 

5.1. Schechter Luminosity Function 

Here we revisit the idea of describing the ob- 
served radio and FIR LFs in terms of Schechter 
functions for two important reasons. First, the 
Schechter function has a simple but sound theoret- 
ical basis as it follows from a theoretical analysis of 
self-similar gravitational condensation in the early 
universe (Press & Schechter 1974). The Schechter 
function has a form 

p{L)dL = p*{L/L*)" exp{-L/L*)d{L/L*) (7) 

where p* and L* are the characteristic density and 
luminosity of the population, and a describes the 
faint-end power-law slope for L <^ L* (Schechter 
1976). Not only is the Press-Schechter formalism 
a well accepted paradigm for cosmological struc- 
ture formation, but the Schechter function also 
describes the observed optical and near-IR LFs 
of field and cluster galaxies well as already men- 
tioned. If star formation activity is an integral 
part of galaxy evolution and reflects the nature 
of the host galaxies, then the Schechter function 
should also offer a good description for the radio 
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and FIR LFs as it does for the optical and near-IR 
LPs. 

Another reason for reconsidering the Schechter 
form is that the observed departure, namely the 
power-law tail at high luminosity end for the FIR 
LF, may have been a predictable consequence of the 
underlying theoretical consideration in the Prcss- 
Schechter formalism. To be specific, the Press- 
Schechter formalism for CDM structure formation 
predicts a mass function of the Schechter form, 
rather than a luminosity function. To make the 
logical connection between the observed LFs and 
the Press- Schechter formalism, one has to invoke 
a constant mass-to- light (M/L) ratio, which is a 
reasonable assumption at the optical and near in- 
frared wavelengths. For starburst galaxies, on the 
other hand, a momentary jump in L/M by 2-3 
orders of magnitudes is expected by definition, 
nearly entirely at the FIR wavelengths. There- 
fore, if the underlying mass distribution is of a 
Schechter form, the resulting FIR (and radio) LF 
should be a linear combination of a Schechter func- 
tion (for field galaxies) and a high luminosity ex- 
cess due to starburst galaxies. Hence, we pose 
a hypothesis that the observed radio and FIR 
LFs are indeed linear sums of two Schechter func- 
tions, one for the field galaxies and second for 
the starburst population. Choosing a Schechter 
form for the starburst population is largely for 
self-consistency. If this population is distinguished 
only by 2-3 orders of magnitude elevation in M/L, 
then a Schechter function should also represent a 
reasonable form for its LF. We examine this hy- 
pothesis by analyzing the derived radio and FIR 
LFs in detail below. 

The characteristic Schechter parameters for the 
observed LFs are derived in two different ways. 
A luminosity function is tlic^ space density of the 
sample galaxies within a luminosity interval of AL 
centered on L. Thus the LF and associated uncer- 
tainty can then be derived "directly" as 

N N 

i=l ''^^ i=l 

where Pm{L) is the measured luminosity density 
within the "magnitude" bin centered on luminos- 
ity L and Vm is the sample volume appropriate 
for each of the N=1809 galaxies (Schmidt 1968; 
Felten 1976). The sample volume Vm was deter- 



mined by the IRAS 60 /im sensitivity limit of 2 Jy 
and the NVSS survey area as: 

/ r^dr dp cos Sd5 = 3 AAD^ . (9) 

Jo Jo Js=-40° 

The sampled volume is further reduced by 18% 
due to the exclusion of the Galactic Plane and the 
omitted areas in the IRAS redshift survey. 

The Schechter parameters can also be derived 
from the observed luminosity distribution, n{L). 
For a randomly chosen sample, the distribution 
function n{L) is given by 

n{L) = p{L)V{L), (10) 

where V{L) is the volume sampled. The distribu- 
tion function can then be re-written as 

n{L)dL = p*V*{L/L*Y+^/^exp{-L/L*)d{L/L*) 

(11) 

where V* is the characteristic volume associated 
with L* . Given these analytical expressions for the 
luminosity distribution and luminosity function, a 

analysis is performed on the observed p{L) and 
n{L) to deduce the characteristic Schechter pa- 
rameters. The Schechter parameters derived from 
p{L) and n{L) are slightly different because the 
two methods are subject to slightly different sys- 
tematic effects. In addition to the formal uncer- 
tainties given by the analysis, the difference be- 
tween the best-fit parameters derived by the two 
different methods offers a quantitative measure of 
the systematic errors in the derivations. 

A potential source of confusion in comparing 
the derived LFs from one study to another is the 
definition of luminosity interval AL used. Both 
"per magnitude" and "per dex" are commonly 
used in the literature. We adopt "per magnitude 
(mag~^)" in our derivations of radio and FIR lu- 
minosity functions. All published luminosity func- 
tions using "dcx^^" intervals (e.g. Saunders ct al. 
1990) are thus corrected by a factor of 2.5 before 
making any comparisons. 

5.2. IRAS 60 pro. Luminosity Function 

The IRAS 60 pm luminosity distribution and 
luminosity function for the IRAS 2 Jy sample 
galaxies are constructed using 18 luminosity bins 
of log 

-^eo/nm (^©) between 7.75 and 12.25 with a 
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width of AlogL ~ 0.25. The computed 60 /im 
luminosity distribution and LF are tabulated in 
Table 4 along with the computed V/Vm, which is 
a measure of the sampling completeness (Schmidt 
1968; Felten 1976). These quantities are also 
shown graphically in Figure 7 along with the best- 
fit Schcchtcr functions (see below). The luminos- 
ity distribution shown in Figure 7a nicely indicates 
that the most common objects in the IRAS 2 Jy 
sample have 60 /im luminosity of about 10^^ Lq, 
typical of L* galaxies in the field. 

The derived 60 fim LF is shown in Figure 7b 
and is compared with previous derivations by 
Soifer et al. (1987) and Saunders et al. (1990). 
The general shape of the LF may be reasonably 
described as a power-law distribution, but two or 
more power-laws are needed to account for the de- 
tailed shape in the entire luminosity range. As dis- 
cussed in the previous section (§ 5.1), a linear sum 
of two Schechter functions may make a good phys- 
ical sense as a description for the radio and FIR lu- 
minosity functions, and indeed the derived analyt- 
ical expressions shown in the solid and dashed lines 
match the the observed 60 /im LF closely. In de- 
riving the best fit Schechter parameters, the faint 
end power-law slope a is assumed to be universal 
because it cannot be derived cleanly for the higher 
luminosity ( "starburst" ) population. The charac- 
teristic luminosities L* derived from the 60 /zm lu- 
minosity distribution (Figure 7a) for the high- and 
low-luminosity populations are (21 ± 2) x 
and (2.2 ± 0.4) x 10^°Lq, respectively, with the 
faint end power-law slope a = —0.81 ± 0.06. 
The best fit parameters can also be derived from 
the observed 60 /im luminosity functions (Fig- 
ure 7b), and they arc (22 ± 2) x 101°Lq and 
(2.3 ± 0.5) X lOi"i0, with a = -0.82 ± 0.05 (see 
Table 5 for the summary and the analysis). 

A break in the FIR LF has previously been sug- 
gested by Smith et al. (1987) and Samidcrs et al. 
(1990) based on the optical morphology and in- 
frared color. These authors had not offered any 
clear physical explanations for such a break seen 
near Leo/jm = lO^^i©, however. The combination 
of an even cleaner break seen in the 1.4 GHz radio 
LF (see below) and the excellent fits reported in 
Table 5 lend strong support for the two Schechter 
function description and the implied physical ex- 
planation in terms of an elevated mass-to-light ra- 
tio for the starburst population. 



Significant differences arc seen among the sev- 
eral 60 yum LFs compared in Figure 7b at the low 
luminosity end. Our 60 /xm LF (filled circles) is 
smoother and has smaller error bars than those de- 
rived by Soifer et al. (1987) and by Saunders et al. 
(1990) because our sample size is 6 times larger. 
Our 60 /im LF dips below the Soifer et al. LF 
(empty squares) and the derived Schechter model 
near ieo^im ^ 10^'^ Lq. The 60 /zm LF derived 
by Saunders et al. (1990; empty circles) fiattens 
even further, beginning around log igo/jm ~ 9.5. 
This fiattening is also seen in their Figure 2, in 
comparison with that of BGS sample (Soifer et al. 
1987), Smith 2 Jy sample (Smith et al. 1987), and 
Strauss & Huchra (1988). One explanation for the 
different faint end slope may be the sample com- 
pleteness as the surveys with particular emphasis 
on the completeness (e.g. Smith et al., Stauss & 
Huchra) tend to derive a steeper faint end slope. 
In contrast, Saunders et al. data probably suf- 
fers the worse completeness problem because their 
luminosity function is derived using a heteroge- 
neous collection of several different surveys. The 
plot of V/Vm (Fig- 7c) shows that our sample is 
distributed fairly uniformly within the sampling 
volume {V/Vm ~ 0.5). 

The derived V/Vm for our sample dips below 
0.5 at log L^Ofim ~ 8.5, albeit with large error- 
bars. This can be interpreted as an indication of 
local density enhancement and may be reflecting 
our special location within the local galaxy group 
and the local super-cluster. Derivation of LFs us- 
ing the 1/Vm method (Eq. 8) has been suggested 
to be susceptible to the large scale inhomogeneity, 
and the steep rise in the derived LF at the faint 
end may be an artifact of this bias. The steep- 
est rise in the faint end slope is associated with 
the LF derived by Soifer et al. which also has the 
smallest sampling volume, in accordance with this 
expectation. Assuming the luminosity function 
has a universal form, Yahil et al. (1991) proposed 
a density-independent determination of the lumi- 
nosity function by introducing a selection function 
which accounts for the distance effect. The Yahil 
et al. derivation of a double power-law form LF 
using an earlier version of the IRAS redshift cata- 
log is shown as a dotted curve in Figure 7b. While 
the derived faint end slope is flatter and shallower 
than the LF we derive, it is also steeper and higher 
than the LF derived by Saunders et al. (also see 
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their Fig. 7). Even such a density-independent 
method faces a fundamental Umitation if the faint 
end slope is constrained by the low luminosity ob- 
jects that are drawn mainly from the local volume 
with enhanced density. Therefore the analysis pre- 
sented here suggests that the faint end behavior of 
the luminosity function is intrinsically poorly con- 
strained. It should be noted, however, that these 
differences between the various LFs at low L do 
not have much effect on the implied star forma- 
tion density discussed below (§6.4) since the star- 
formation rate is proportional to the LF multiplied 
by L, a function peaking at higher luminosity. 

5.3. 1.4 GHz Radio Luminosity Function 

Because our sample is selected by the IRAS 60 
fim band flux density, deriving the 1.4 GHz radio 
luminosity function for the same sample is not as 
straight forward as in the 60 /xm LF case. In par- 
ticular, the 1.4 GHz radio LF cannot be derived 
from the luminosity distribution (Eq. 11) because 
the volume sampled does not relate directly to the 
1.4 GHz luminosity. Instead, the 1.4 GHz radio 
luminosity function is derived using the sampling 
volume corresponding to the IRAS 60 /xm luminos- 
ity of each sample galaxy. The logic behind this 
calculation is that each sample galaxy contributes 
to the local number density at their respective lu- 
minosity bins in exactly the same manner for both 
the 60 /zm and 1.4 GHz LF, and sample selection 
is reflected identically by the sample volume de- 
termined by the IRAS 60 flux limit. 

The computed radio luminosity distribution 
and LF are tabulated in Table 6 along with the 
computed V/Vm. These quantities are also shown 
graphically in Figure 8 along with the best-flt 
Schechter functions (see Table 7). The radio lu- 
minosity distribution, LF, and the V/Vm plots 
are qualitatively similar to those of the IRAS 
60 jum functions in Figure 7, mainly because of 
the linear radio-FIR correlation (see §3 & §6.3). 
The derived 1.4 GHz radio LF, shown in Fig- 
ure 8b, is flat at low luminosity end and falls off 
steeply beyond 10^^ W Hz~^. A Schechter func- 
tion (solid line) appears to be a good description 
for all low luminosity objects of Li^ghz < 10^^ 
W Hz~^, with the best fit characteristic Schechter 
parameters of L* = (2.1 ± 0.3) x 10^2 W Hz"!, 
p* = (3.2 ± 0.2) X 10-* Mpc-3 mag-\ and 
a = —0.63 ±0.05. The observed radio LF is signif- 



icantly flatter than the best fit Schechter function 
below LiAGHz < 10^°-* W Hz-i. In addition to 
the possible effects of sample completeness and 
local density enhancement, the systematic devia- 
tion from the linear radio-FIR relation may further 
contribute to this flattening. 

The presence of the luminous second popula- 
tion is more obvious in Figure 8b at Li^ghz > 
10^3 w Hz"^. A fuU analysis is not war- 
ranted for this population because of the small 
numbers of bins represented and the large ob- 
served scatter. Characteristic luminosity and den- 
sity of L* = (1.4 ± 0.2) X 10^3 W Hz-i and 
p* = (8.3 ± 0.8) X 10"*^ Mpc"^ mag-i are derived 
assuming the same faint-end power law slope of 
a = —0.63, but these mmibers should be taken 
only as illustrative values. The flux cutoff of 2 
Jy at 60 ^m and the resulting volume sampled 
are probably not sufficient to fairly represent these 
rarest objects with high luminosity. 

The comparison of our radio LF (filled circles) 
with that of the flux hmited sample of UGC galax- 
ies by Condon (1989; open circles) shows that 
the low radio luminosity galaxies in our IRAS se- 
lected sample are nearly identical to the "star- 
burst" galaxies (late type galaxies in the field, 
as opposed to "monsters") in the UGC catalog 
(see Figure 8b). The small excess associated with 
the IRAS 2 Jy sample is entirely accounted by 
the cluster galaxy population. When all galaxies 
within the 3 Mpc projected diameter of the known 
clusters within 100 Mpc distance are removed, the 
resulting radio LF is indistinguishable from the 
UGC sample LF. The significance of this result is 
that there is a complete overlap between the opti- 
cally selected late type galaxies and IRAS selected 
galaxies in the field. In other words, the measured 
radio luminosity may be used to infer star-forming 
activity for late-type galaxies in the field with the 
same level of confidence as using their FIR lumi- 
nosity. Estimating the star-formation rate for an 
individual galaxy using its radio luminosity may 
still be uncertain by the magnitude of the observed 
scatter in Figure 6, but such an analysis for an 
ensemble of galaxies should be significantly more 
reliable (see §6.5). 

The comparison with the "monsters" in the 
UGC sample (open squares), which has a flat lu- 
minosity function, is more complicated. The 1.4 
GHz radio LF for the "monsters" derived by Con- 
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don is comparable to that of our IRAS selected 
galaxy sample for L\aghz between 10^'^ and 10^'' 
W Hz~^, but the difference quickly grows to more 
than an order of magnitude at Liuqhz <^ lO^** 
W Hz~^. This disagreement at the highest lu- 
minosity is highly statistically significant, and it 
must be explained by the incompleteness in the 
IRAS 2 Jy sample. The LF of galaxies is a prop- 
erty of galaxies and not of the sample; the 
weighting is supposed to correct for sample selec- 
tion effects. With a sufficiently large sample of 
IRAS galaxies, the LF should equal the LF of op- 
tical/radio selected UGC galaxies, including both 
starbursts and monsters. The main reason that 
our IRAS total LF is lower than the UGC mon- 
ster LF must be that NO "monsters" have been 
detected by IRAS in the highest luminosity bins. 
If even one or two were detected, they would have 
small V„i values and hence large weights, bring- 
ing up the IRAS LF to agree with the UGC LF. 
Therefore, the second, more luminous population 
in the IRAS 2 Jy sample cannot be simply iden- 
tified with the "monsters" in the UGC sample. 
Because these high luminosity objects also largely 
follow the same radio-FIR correlation as the low 
luminosity objects (see Fig. 5), the observation can 
be explained more naturally with our hypothesis 
that the high-luminosity excess represents a phys- 
ically distinct second population, whose large lu- 
minosity is the manifestation of a burst of tidally 
or merger induced activities (sec § 5.1). Possible 
contribution by luminous AGNs are discussed in 
greater detail below. 

6. Discussion 

6.1. High Luminosity Objects and Fre- 
quency of AGNs 

Infrared-luminous galaxies with Li?/fl > lO"^""^ Lq 
(equivalently, Li^ghz > 10^'^ W Hz~^) are spe- 
cial objects because few normal disk galaxies in the 
local universe have such large luminosity. A com- 
plete list of 161 luminous infrared galaxies with 
Leofi > W^^'^Lq is given in Table 1. In addition to 
many well-known ultraluminous infrared galaxies 
(e.g. Arp 220), Seyferts, starburst, and HII galax- 
ies are included in the list, and morphological 
descriptions of "interacting", "galaxy pairs", or 
"mergers" are common. More importantly, they 



may represent a distinct separate population from 
the late type field galaxies in the derived radio 
and FIR LFs (see §5). 

The infrared color-color diagrams shown in Fig- 
ure 9 provide a useful tool for probing the mean 
radiation field and dust composition. In analyz- 
ing the IRAS colors of field galaxies, Helou (1986) 
found a linear trend of increasing dust tempera- 
ture with increasing luminosity as shown by the 
dashed parallel tracks in Figure 9. The high- 
luminosity objects (filled large circles) as a group 
show a systematically larger 60 ^m to 100 fim flux 
ratio, which is an indication of warmer mean dust 
temperature {Seonm/Siooum ~ 0.8, Tdust ~ 45 K). 
They also cluster near the upper left corner of 
Figure 9, which indicates that the dust in these 
galaxies is exposed to hundreds of times stronger 
radiation field than that of the solar neighbor- 
hood. This clustering is broadly consistent with 
the intense starburst interpretation for these FIR- 
luminous systems. Even though their FIR lumi- 
nosity is on average 10 times smaller, the radio- 
excess objects (radio AGNs, shown as diamonds) 
show a similar IRAS color distribution as infrared 
luminous galaxies. An AGN produces its charac- 
teristic mid-IR enhancement by heating a small 
amount of dust in its immediate surroundings to 
a very high temperature. The observed similarity 
in the IR color between luminous starbursts and 
AGNs might indicate that the warm dust heated 
by stars may also be concentrated in a relatively 
small volume. 

While an unambiguous distinction is not possi- 
ble, a systematic separation among the late type 
field galaxies, infrared luminous galaxies, and ra- 
dio AGNs becomes more apparent in the far-IR 
versus mid-IR color plot as shown in Figure 10. 
The far-IR color of the FIR-luminous sample is 
clearly and systematically warmer than that of the 
late type field galaxies (small dots). The median 
Seoum/ SiQOfim ratio for the FIR-luminous galaxies 
is about 0.8 (45 K for /? = 1; see Helou et al. 1988) 
while it is about 0.5 (36 K for /3 = 1) for the field 
galaxies. The radio-excess objects show a similar 
range of the far-IR colors as the infrared luminous 
objects. In addition, they also show a systematic 
offset towards the warmer mid-IR color by about 
0.5 in dex, separating themselves from the infrared 
luminous galaxies. Presence of hot dust in the 
circumnuclear region surrounding an AGN offers 
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a plausible explanation for this mid-IR enhance- 
ment. Mid-IR excess with S25nm/ S^oum > 0.18 
(right of the thin vertical line) has been empirically 
proposed as an indicator of an infrared AGN by 
de Grijp et al. (1985). Indeed the radio-excess ob- 
jects in our sample show the largest S2hixm/ S^oum 
ratios among all IRAS selected galaxies (see Fig- 
ure 10). 

A total of 20 out of the 161 FIR-luminous ob- 
jects in our IRAS 2 Jy sample have S25fj.m/ Seofim > 
0.18. This fraction (20/161=12%) can be inter- 
preted as an upper limit to the fraction of FIR- 
luminous galaxies that may be powered primarily 
by an AGN. Genzel et al. (1998) found that an 
energetically dominant active nucleus is present 
in 20 to 30% of all FIR-ultraluminous galaxies 
based on the analysis of their mid-IR diagnostic 
diagram. This broad agreement is not surprising 
since both methods utilize the enhanced mid-IR 
emission as an indicator. The actual fraction of 
FIR-luminous galaxies powered by AGNs may be 
smaller still as discussed below. The total fraction 
of radio-excess objects among the mid-IR excess 
objects is only 10% (2 out of 20). Perhaps for- 
tuitously, this fraction is similar to the frequency 
of radio-loud objects among optically identified 
QSOs (Kellermann et al. 1989; Hooper et al. 
1995). 

The proposed separation of AGNs from star 
forming galaxies using the »S'25/^m/*S*6o^m flux ra- 
tio is not entirely supported by Figure 10, how- 
ever. Both the field galaxies and FIR-luminous 
galaxies show a broad general trend towards a 
higher S25fmi/ Seoiim flux ratio with increasing 
SeOfim/ SiQOfim flux ratio. This trend is seen more 
clearly when each class of objects are considered 
separately, and this follows the expected behav- 
ior for an aging and evolving starburst. To quan- 
tify this trend, the color evolution track for a 
starburst galaxy computed by Efstathiou et al. 
(2000) is shown as a thick curve in Figure 10. 
In this model, a starburst at its onset (top right 
end of the track) displays significant mid-IR emis- 
sion contributed by the hot /warm dust directly 
exposed to the intense radiation field of the young 
massive stars. The predicted S'25^m/'5'60/xm ratio 
is comparable to the largest values observed for 
the FIR-luminous galaxies. As the starburst ages 
and the hot /warm dust cools, the starburst trav- 
els down along the track as shown in Fig. 10. It 



is particularly interesting that this color evolu- 
tion track closely hugs the proposed AGN sepa- 
ration line at S25fim/Sc,otJ.m = 0.18 and outlines 
the right side boundary for the FIR-luminous ob- 
jects in this plot. The treatment of the super- 
nova ejecta in this model is such that the evolu- 
tion track shown is essentially the right boundary 
for a family of possible evolution tracks, in agree- 
ment with the data shown here - see Efstathiou et 
al. (2000) for a detailed description of the model. 
The youngest starbursts (< 20 Myr) are predicted 
to appear well to the right of the AGN separa- 
tion line, and this prediction is further supported 
by the fact that young starburst galaxies hosting 
massive Wolf-Rayet clusters such as NGC 1614, 
NGC 1741, and NGC 5253 show Scyfcrt-like mid- 
IR excess with S2bixm/ S^oum flux ratios of 0.22, 
0.15, and 0.39, respectively, while they show no 
evidence for hosting an AGN. Therefore some cau- 
tion should be taken for interpreting mid-IR excess 
as an AGN indicator, especially for galaxies with 
warm far-IR color. 

The fact that an overwhelming majority of 
high-luminosity objects in our IRAS selected sam- 
ple do not show mid-IR excess and have the same 
infrared color characteristics as the less luminous 
normal star forming galaxies would argue that 
most of them are powered by an intense star- 
burst instead. The derived radio LF for the high- 
luminosity objects in the IRAS 2 Jy sample be- 
ing not consistent with that of the radio AGNs or 
"monsters" in the UGC sample of Condon (1989) 
also indicate that the majority of them may not be 
powered by an AGN. Further, about 30% (7/23) of 
the radio-excess objects appear to the left of the 
mid-IR excess line in Figure 10, indicating that 
some of the radio AGNs do not contribute signifi- 
cantly to their host's infrared luminosity. 

6.2. Nature of FIR-Excess Galaxies 

There arc nine objects in the IRAS 2 Jy sam- 
ple with unusually large FIR/radio flux ratios with 
g > 3. Three objects have ^eo^m/S'iooMm > 
1 {Td = 60 ^ 80 K), and the dust temper- 
atures in these galaxies are among the highest 
in the sample. The strong temperature depen- 
dence in dust cmissivity implies that a dramatic 
increase in FIR luminosity can result from only 
a slight increase in dust temperature. Similarly 
large SQOnm/Siooum ratios are also found in in- 
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frared quasars FSC 09105+4108 (Hutchings & 
Neff 1988), FSC 13349+2438 (Beichman et al. 
1986), FSC 15307+3252 (Cutri et al. 1994) and 
other luminous AGN host galaxies (see Yun & 
Scoville 1998). Ironically, broad emission lines 
indicative of a luminous AGN have not been de- 
tected by optical or infrared spectroscopy among 
the FIR-excess galaxies observed to date. Ev- 
idence for large visual extinction and/or large 
gas and dust contents is present in the four best 
studied galaxies (NGC 4418, IRAS 08572+3915, 
IRAS 20087-0308, IRAS 22491-1808). 

Using the high resolution imaging of 8.4 GHz 
radio emission in the 40 ultraluminous galaxies in 
the IRAS Bright Galaxy Sample, Condon et al. 
(1991b) have shown that starburst regions in some 
of these galaxies are compact and dense enough 
to be opaque even at mid-IR and long radio wave- 
lengths. Such compact and intense starbursts can 
indeed account for the large Seo^m/ Siooum ratios 
and high extinction at the same time. These in- 
frared excess objects may actually represent a spe- 
cial phase in the starburst evolution. 

The existing data, however, cannot uniquely 
distinguish whether these FIR-excess objects are 
powered by a dust enshrouded AGN or by a com- 
pact starburst. If a dense environment and associ- 
ated large opacity is responsible for the observed 
infrared excess, then one cannot possibly probe 
the energy source directly since the central active 
region must be opaque - even for hard X-ray pho- 
tons if the line of sight column density exceeds 10^'' 
cm~^. If free- free absorption is responsible for the 
low 1.4 GHz flux, then dense ionized gas must fill 
most of the 1.4 GHz source volume (> 50 — 100 pc 
in size). Since the UV opacity is so high, an ion- 
izing photon cannot travel very far before being 
absorbed by dust. Thus an extended, distributed 
source (a starburst) is favored over a single point 
source (an AGN) to account for the observed ion- 
ization. 

What is more certain is that such a compact 
starburst phase or a dust-enshrouded AGN phase 
must be brief, perhaps only a few percent of the IR 
bright phase, because only nine out of 1809 sources 
show such infrared excess.^ Converting significant 

'Although this cutoff is somewhat arbitrary, they represent 
the most extreme cases among the galaxies responsible for 
the symmetric dispersion in Figure 6, and such infrared- 



fraction of the AGN luminosity into thermal dust 
emission requires a geometry for thoroughly en- 
shrouding the AGN with dust. The paucity of 
FIR-excess objects in the IRAS 2 Jy sample thus 
argues strongly that FIR-luminous galaxies pow- 
ered by luminous AGNs account for a very small 
fraction. The AGN-powered FIR-luminous galax- 
ies may be in the radio-excess phase instead, but 
we have already established in § 4.2 that the radio- 
excess fraction is less than a few per cent, indepen- 
dent of luminosity. 

6.3. Evidence for Non-linearity in the 
Radio-FIR Correlation 

As a further test of the radio-FIR correlation, 
the IRAS 60 fiia LF is plotted on top of the 1.4 
GHz radio LF in Figure 11 after shifting the 60 fim 
LF along the x-axis using the linear radio-FIR re- 
lation in Eq. 4. The agreement is extremely good 
for the two LFs over a wide range of luminosity. A 
total of 1467 galaxies (81%) belong to the luminos- 
ity range between 10^^ and 10^^ W Hz^^ where the 
two LFs are essentially identical. The luminosity 
range over which the two LFs agree within a fac- 
tor of two includes over 98% of the total sample. 
Thus, the similarity between the two LFs stems 
from the fact that the overwhelming majority of 
galaxies in our IRAS selected sample follow the 
same linear radio-FIR correlation. 

Evidence for diminished radio emission among 

low-luminosity galaxies was discussed earlier (see 
§3.1 and in Fig. 6), but its collective impact on the 
radio and FIR LF appears to be minimal. This 
non-linearity among the low-luminosity galaxies 
should result in flattening of the radio LF with 
respect to the 60 /xm LF. Indeed a hint of such 
an offset is seen in Figure 11, but the difference is 
marginal. A more substantial difference is seen at 
the high luminosity end, where a dramatic increase 
in the radio LF over the 60 /im LF is seen near 

LlAGHz ~ 1023-5 W Hz-1 (L60/.r» ~ IO^^'^Lq). 

While this may be interpreted as a sign of in- 
creased AGN activity associated with the high lu- 
minosity population, the falloff in the LF by an 
order of magnitude near Liaghz ~ 10^* W Hz~^ 
is not consistent with such a supposition (see §5.2). 
The frequency of radio-excess objects was found to 

excess objects are extremely rare. 
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be independent of luminosity in §4.2, and the con- 
stant AGN fraction should appear as a small and 
constant offset between the two LFs. A plausible 
explanation for the large scatter may be the small 
number statistics in the total numbers of objects 
and the AGN frequency in these bins. 

6.4. FIR Luminosity Density and Local 
Star Formation Rate 

The derivation of the integrated IRAS 60 /xm 
luminosity density for the local volume is shown in 

Figure 12. The derived 60 /im luminosity density 
integrated between Leo/im = 

lO^L© and lO^^L© 
for the IRAS 2 Jy sample is 2.6 ± 0.2 x lO'^L© 
Mpc^"*. This estimate lies midway between the 
values derived from the 60 /um LFs of Soifer et 
al. (1987) and "S17" of Saunders et al. (1990). 
These differences arise mainly from the different 
behaviors of these LFs at low luminosity end (see 
Fig. 7b and §5.3), and the observed scatter can be 
interpreted as a measure of uncertainty resulting 
from the sample selection and clustering of galax- 
ies. In all three cases, the integrated luminosity 
density reaches the asymptotic peak value near 
Leoum = lO^^I'o, and the infrared luminous galax- 
ies (LpiR > 10^^ Lq) contribute less than 10% 
to the FIR luminosity density in the local volume 
(z < 0.15). 

Using Eq. 3 and assuming a median IRAS 60 
lim to 100 /im flux ratio of 0.5 (see Figure 10), we 
estimate a local FIR luminosity density of (4.8 ± 
0.5) X lO'^Z/0 Mpc~^. This is midway between 
the values derived by Soifer et al. (~ 9 x IO'^Lq 
Mpc~3) and Saunders et al. (4.2 x lO^L© Mpc^^). 
The derived FIR luminosity density can then be 
used to derive the extinction- free star formation 
rate (SFR) for the local volume assuming the ob- 
served FIR luminosity is proportional to the level 
of star forming activity. 

A wide range of conversion relations between 
SFR and LpiR are found in the literature, based 
on either a starburst model or an observational 
analysis (e.g., Hunter et al. 1986, Meurer et al. 
1997, Kennicutt 1998ab). Using the continuous 
burst model of 10-100 Myr duration by Leitherer 
& Heckman (1995) and Salpeter IMF with mass 
Umits of 0.1 and 100 Mq, Kennicutt (1998ab) de- 
rives a conversion relation appropriate for star- 



bursts as 

SFR (M^yr-') = 1.7 x lO-'^'Ls-woOf^miLe) 

(13) 

with about 30% uncertainty. The ratio between 
Ls-wooi^m and LpiR {= L40-120mto) for the 60 
IRAS BGS galaxies studied by Sanders et al. 
(1991) ranges between 1.15-2.0 with a median 
value of about 1.30. Based on actual measure- 
ments at both long and short wavelengths, Meurer 
et al. (1999) find the ratio is larger than 1.4 while 
Calzetti et al. (2000) find an average value of 
1.75 for the five galaxies they examined in de- 
tail. Adopting a conservative ratio of 1.5 and us- 
ing Eq. 13, we derive the star formation rate of 
0.012 ± 0.004 Mq yr-i Mpc'^ for the local vol- 
ume {z ^ 0.15). Since more than 90% of the local 
FIR luminosity density is contributed by normal, 
late-type field galaxies rather than by the star- 
burst population, a conversion relation derived for 
late-type field galaxies may be more appropriate 
for the IRAS 2 Jy sample. The local star forma- 
tion density derived using the conversion relation 
derived for the field galaxies by Buat & Xu (1996) 
is slightly larger, 0.015 ± 0.005 Mq yr"! Mpc'^. 
Contribution from cold cirrus and smaller optical 
depth in late type galaxies add some uncertainty 
to the inference of star-formation rate from FIR 
luminosity (see Kennicutt 1998a for a detailed dis- 
cussion). 

The two different methods produce a consistent 
estimate for the local star formation density. In 
comparison, the the local SFR derived from the 
Ha luminosity density measured by Gallego et al. 
(1995) using the conversion relation given by Ken- 
nicutt (1998a) is 0.014 ± 0.007 Mq yr'^ Mpc'^. 
An important consequence is that the local SFR 
derived from the FIR luminosity density provides 
an independent estimate for 2; = and lends a 
further support for the z = normalization for 
the cosmic star-formation history models such as 
those by Blain et al. (1999; dashed line) and by 
Tan et al. (1999; solid line) as well as support for 
the larger extinction correction for all UV derived 
values as proposed by Steidel et al. (1999) and 
others (see Figure 13). 

The SFR derived from the FIR luminosity den- 
sity may be an overestimate if the observed FIR in- 
cludes a significant contribution from dust heated 
by an AGN or by interstellar radiation field. The 
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frequency of the radio-excess objects in our sample 
is only about 1%, independent of luminosity. And 
the overall frequency of AGNs (including radio- 
quiet) may be as high as 10% if only one out 
of ten AGNs produces significant radio emission. 
Some AGN contribution to the total FIR lumi- 
nosity density is undoubtedly present, but it is 
nevertheless not likely to be dominant (see §6.1 & 
§6.2). The contribution from the dust heated by 
interstellar radiation field has been suggested pre- 
viously to account for up to 50 to 70% of the FIR 
luminosity (see Lonsdale-Persson & Helou 1987, 
Sauvage & Thuan 1992, Xu & Helou 1996). This 
claim does not seem entirely consistent with our 
finding that an overwhelming majority of both the 
field and starburst galaxies obey the same linear 
radio-FIR relation (see §3 and §6.2). A possible 
explanation may be that the cirrus is heated by 
non-ionizing UV photons from relatively massive 
stars, so that even the cirrus tracks the population 
of stars producing the radio emission (Xu 1990). 

6.5. 1.4 GHz Radio Luminosity and Star 

Formation Rate 

The excellent linear correlation seen between 
the FIR and radio luminosity among the IRAS 2 
Jy sample offers a possibility of deriving the star 
formation rate directly from the measured radio 
luminosity. Because of the significant scatter in 
the radio-FIR relation and possible non-linear ef- 
fects discussed in §3, determining the star forma- 
tion rate for any particular galaxy using its radio 
luminosity alone may be risky. The effect of the 
scatter is substantially reduced, however, when 
the radio LF is integrated over luminosity in de- 
riving the radio luminosity density. Further, the 
star formation density for a given volume can be 
determined with a greater accuracy if the integra- 
tion is limited over the range of radio luminosity 
where the linear radio-FIR relation is secure. 

Deriving the star-formation density from the 
radio continuum luminosity density for the general 
galaxy population also has to account for the pos- 
sible contribution from radio AGNs. As shown in 
Figure 14, the 1.4 GHz radio LF is dominated by 
the normal late type galaxies for Lxaghz < 10^^ 
W Hz~^ while the contribution by the "monsters" 
becomes significant only at Liaghz > 10^** W 
Hz~^. This bimodality is also clearly seen in the 



histogram of 1.4 GHz radio luminosity distribu- 
tion for Condon's UGC sample shown in Figure 15. 
Since the net AGN contribution is small and quan- 
titatively well imderstood for Li ^qhz < 10^^ W 
Hz~^, the star formation density for the local vol- 
ume can be derived with some confidence from the 
radio luminosity density integrated up to this lu- 
minosity cutoff. 

The integrated 1.4 GHz luminosity density de- 
rived from the IRAS 2 Jy sample for Li ^gHz < 
10^4 W Hz-i is 2.6 X IQi^ W Hz"! Mpc'^ (see 
Figure 14). The empirical conversion relation that 
links this 1.4 GHz luminosity density to the local 
star formation density of 0.015 ± 0.005 Mq yr~^ 
is then 

SFR (MQyr-^) = 5.9±1.8xl0-^2^i iGHz{W Hz-^). 

(14) 

As discussed in §6.4, the main micertainty in 
this relation comes from the estimate of the local 
star formation density. If the radio LF at a given 
epoch can be derived with sufficient accuracy, then 
the contribution by star forming galaxies and radio 
AGNs can be distinguished and the star formation 
density for the epoch can be derived using Eq. 14. 
Whether the dividing line between the starburst 
and AGN population near Liaghz ^ 10^^-^ W 
Hz""'^ is constant or evolved over time is unknown 
at the moment, and it needs to be investigated by 
future studies. Ultimately, the cosmic star forma- 
tion history free of any extinction effects may be 
derived using the radio emission from star form- 
ing galaxies and by tracking the evolution in their 
radio luminosity function. 

7. Conclusions 

By cross-correlating the IRAS Redshift Survey 
catalog and the NRAO VLA Sky Survey database, 
we have assembled a flux-limited complete sam- 
ple of 1809 IRAS detected galaxies. This IRAS 
2 Jy sample is 6 times larger in size and 5 times 
deeper in redshift coverage {z < 0.15) than pre- 
vious studies of radio and infrared properties of 
galaxies. In addition to obtaining the accurate 
positions of these IRAS sources, we have exam- 
ined the well known radio-FIR correlation in de- 
tail. The principal conclusions from the analysis 
of the data are the following: 

1. The radio-FIR correlation for the FIR- 
selected galaxy sample is well described by a linear 
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relation, and over 98% of the sample galaxies fol- 
low this linear radio-FIR correlation. The scatter 
in the linear relation is about 0.26 in dex, domi- 
nated by small systematic deviations involving a 
small number of galaxies at high and low luminos- 
ity ends. 

2. A total of 23 "radio-excess" objects (radio 
AGNs) are found in our sample including three 
"radio-loud" objects PKS 1345-1-12, 3C 273, and 
NGC 1275. The frequency of radio-excess ob- 
jects in our IRAS 2 Jy sample is 23/1809 = 1.3%, 
and radio-excess objects are rare among the FIR- 
selected sample of normal and starburst galaxies. 
There are only 3 radio-excess objects among the 
161 FIR-luminous (ieo/im > lO^^'^L©) galaxies in 
the sample (3/161 = 1.9%). Therefore the fre- 
quency of radio AGNs appears to be independent 
of FIR luminosity. 

3. The derived 1.4 GHz radio luminosity func- 
tion and 60 /im luminosity function for the IRAS 
2 Jy sample are both reasonably modeled as linear 
sums of two Schechter functions, one representing 
late type field galaxies and the second represent- 
ing a starburst population. The presence of the 
second, high luminosity component in these LFs 
is a predictable outcome of the elevated M/L as- 
sociated with the starburst population, and the 
Press-Schechter formulation of the structure for- 
mation can now be naturally extended to the in- 
terpretation of the radio and infrared LFs as well. 

4. The radio-excess objects in our sample show 
evidence for the same highly enhanced radiation 
field as the infrared luminous galaxies in the in- 
frared color-color diagrams. A systematic segrega- 
tion of field galaxies, starburst galaxies, and radio 
AGNs is found in the mid-IR versus far-IR color- 
color diagram (Fig. 10), where the radio AGNs 
show the largest observed S2r->tim/ SaQfj^m ratios. 
The total fraction of radio-excess objects among 
the galax;ies with mid-IR enhancement is only 
about 10% (2/20) - similar to the frequency of 
radio-loud objects among the QSOs. The absence 
of mid-IR enhancement and radio-excess among 
the large majority of the IR luminous galaxies sug- 
gests that star formation is the dominant energy 
source in most cases. 

5. A small number of infrared-excess objects 
are also identified by their departure from the 
radio-FIR relation (9 objects with q > 3). Both 
dust-enshrouded AGNs and compact starbursts 



can explain their observed properties, and distin- 
guishing the two may be difficult. The observed 
spatial extents of 50-100 pc favors the starburst 
explanation. In either case, this phase must be 
brief, perhaps lasting only a few percent of the IR 
luminous phase since these infrared-excess objects 
are extremely rare. 

6. The integrated 60 jim luminosity density 
between Leo^m = lO^L© and IO^^Lq is 2.6 ± 
0.2 X W'^Lq Mpc--\ The FIR-luminous galax- 
ies {LpiR > lO^^L©) contribute less than 10% of 
the total. The derived local FIR luminosity den- 
sity is (4.8 ± 0.5) X IO^Lq Mpc'^ The inferred 
extinction-free star-formation density for the local 
volume is 0.015 ± 0.005 Mq yr'^ Mpc'^. 

7. The observed linear correlation between the 
FIR and radio luminosities in the IRAS 2 Jy sam- 
ple offers a potentially useful method of deriving 
the star-formation rate using the measured radio 
luminosity density. A good separation found be- 
tween star-forming galaxies and AGNs in the de- 
rived radio LF offers a promising possibility of 
tracing the evolution of the two populations sep- 
arately by deriving their radio LFs at different 
epochs. 

By examining the radio properties of a large 
sample of infrared selected galaxies, we have ob- 
tained a much better understanding of the ex- 
tent of the radio-FIR correlation. We have also 
obtained a good quantitative understanding of 
the frequency and energetic importance of active 
galactic nuclei among the infrared selected galax- 
ies and outlined the practical limits on using the 
radio continuum as a direct probe of star forming 
activity. This study can now serve as the basis 
for conducting deep radio, IR/submm, and opti- 
cal surveys aimed at mapping the extinction-free 
cosmic-star formation history and the evolution of 
radio AGNs. 
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Fig. 1. — Equal area sky projection of the IRAS 
2 Jy Sample. Large gap to the right and left of 
the center trace the Galactic Plane, and a narrow 
empty strip near the center is the well known gap 
in the IRAS survey. Otherwise, the sky coverage 
is fairly uniform. 



Fig. 2. — Plot of 1.4 GHz radio luminosity ver- 
sus redshift for the IRAS 2 Jy sample. The lower 
bound is defined by the 2 Jy flux limit of the IRAS 
60 fim band. The upper bound is dictated by the 
radio luminosity function. The three galaxies with 
LiAGHz > 10^^ W Hz~^ are identified as "radio- 
loud" objects. 
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Fig. 3. — An example of confusing source identi- 
fication using the NVSS catalog. The IRAS PSC 
position for IRAS 12173—3541 and its error ellipse 
are centered between the two interacting galaxies 
while the dominant radio source is identified with 
the northern galaxy by the NVSS (marked with a 
white cross). 
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Fig. 4. — A histogram of the offsets in arcsec be- 
tween the IRAS PSC and the NVSS catalog po- 
sitions. The distribution is not Gaussian because 
the IRAS position error ellipse is typically highly 
elongated. The position offset is less than 10" in 
over 2/3 of the cases, which is consistent with the 
expectation from the IRAS position uncertainty. 
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Fig. 5. — (a) Plot of 1.4 GHz radio luminosity 
versus IRAS 60 fj,m luminosity. The solid line cor- 
responds to a linear relation with a constant offset, 
(b) Plot of 1.4 GHz and IRAS 60 /im flux density 
for the IRAS 2 Jy sample. The solid line corre- 
sponds to the same linear relation shown in (a). 
Gray filled circles identify the radio-excess objects 
(see §4.2). The remaining ^1750 objects (out of 
1809) lie very close to the linear relation, and the 
rms scatter of the data is less than 0.26 in dex. 



Fig. 6. — Distribution of q values plotted as a 
function of IRAS 60 /im luminosity. The solid line 
marks the average value of g = 2.34 while the 
dotted lines delineate the "radio-excess" (below) 
and "infrared-excess" (above) objects, delineated 
for having 5 times larger radio and infrared flux 
density than the expected values from the linear 
radio-FIR relation, respectively. 



24 













(a) : 










\ i 




100 








i 














Cou 






A 








10 


r 








o 






i / 








1 






. < 1 . . < . 


.1 ■ 



8 9 10 11 IS 



e 



-2 




(b) : 


-3 






-4 






-5 




^ - 


-6 










\ \ 


-7 


□ Soifer et al. 


(1987) \ \R . - 




o Saunders et 


al, (1990)'^ ff; 


-8 




\, , A* 



B 9 10 11 13 



log Lb(i„„(Lo) 



(c) - 




- 



a g 10 11 12 
log ho^JK) 

Fig. 7. — (a) 60 fim luminosity distribution for 
the IRAS 2 Jy sample. The data points are num- 
ber of counts in each luminosity bin, and the error 
bars reflect the Poisson statistics, N^^^ (see Ta- 
ble 4). The two solid lines indicate the best-fit 
curves of the Schechter form listed in Table 5. (b) 
The 60/im infrared luminosity function. Filled cir- 
cles represent the LF for the IRAS 2 Jy sample. 
The open squares and open circles represent the 
60 fim LFs derived from the IRAS BGS sample 
by Soifer et al. (1987) and an ensemble of data 
("S17") analyzed by Saunders et al. (1990), re- 
spectively. The long dashed lines correspond to 
the two Schechter fits for the high and low lumi- 
nosity objects listed in Table 5. The solid line 
is the resulting linear sum of the two Schechter 
functions. The dotted curve represents the dou- 
ble power-law LF derived from an earlier version 
of the IRAS redshift survey sample by Yahil et al. 
(1991) using a density-independent determination 
method, (c) The V/V^ plot for the IRAS eO^im 
luminosity distribution. The V/Vm analysis sug- 
gests that the highest and the lowest luminosity 
bins are not well sampled and may suffer a com- 
pleteness problem. 
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Fig. 8. — (a) 1.4 GHz radio luminosity distribu- 
tion for the IRAS 2 Jy sample. The data points 
are number counts in each luminosity bin, and the 
error bars reflect the Poisson statistics, N^^'^ (see 
Table 6). The two solid lines indicate the best- 
fit curves of the Schechter form listed in Table 7. 
(b) The 1.4 GHz radio luminosity function. Filled 
circles represent the IRAS 2 Jy sample while the 
open circles and squares represent the 1.4 GHz LFs 
for the "starburst" and "monsters" from Condon 
(1989), respectively. The solid line is the best fit 
Schechter function for the low luminosity objects 
only (see Table 7). (c) The V/Vm plot for the 1.4 
GHz radio luminosity distribution. The error bars 
are statistical and may not be accurate for bins 
with small number of counts. The V/Vm analysis 
suggests that the highest and the lowest luminos- 
ity bins are not well sampled and may suffer a 
completeness problem. 
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Fig. 9. — IRAS mid-IR to far-IR color-color di- 
agram for our IRAS 2 Jy Sample. Large filled 
circles correspond to infrared luminous galaxies 
with ieoM > 10"-^Lo (see Table 1). The low lu- 
minosity, late type galaxies in the field are plot- 
ted in small dots. The large diamonds identify 
the "radio-excess" objects (radio AGNs; Table 2). 
The twin broken parallel lines trace the area occu- 
pied by "normal galaxies" as identified by Helou 
(1986), and the curved solid line is the theoreti- 
cal trajectory of mixed composition dust under in- 
creasing mean radiation field calculated by Desert 
(1986), from solar neighborhood intensity ("cir- 
rus" , lower right) to several hundred times larger 
radiation field (upper left). 



Fig. 10. — IRAS far-IR versus mid-IR color-color 
plot for the IRAS 2 Jy sample. The S'eoAim/'S'iooMm 
ratio is translated into dust temperature assuming 
emissivity index (3=1 (see Helou et al. 1988). All 
the symbols are the same as in Figure 9. The thin 
vertical line marks the proposed division for the in- 
frared Seyferts (S'25Aim/S'60/im > 0.18) by de Grijp 
et al. (1985). The thick solid line represents the 
color evolution track for a model starburst galaxy 
from Efstathiou et al. (2000), starting at zero age 
on the top right corner to 72 Myr old starburst at 
the end (see their Fig. 2 for more detail). 
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Fig. 11.— The 60 /im FIR LF (in filled squares) 
is translated using the linear radio-FIR relation in 
Eq. 4 and compared directly with the 1.4 GHz ra- 
dio LF (in filled circles). The two Schechter func- 
tions for the high and low luminosity 60 /im FIR 
LFs are shown in long dashed lines and the sum 
of the two are shown in thick solid line as in Fig- 
ure 7b. 



Fig. 12. — A plot of the cumulative IRAS 60 /im 
luminosity density. The filled circles mark the 
IRAS 2 Jy sample while the filled squares and 
crosses represent the IRAS BGS sample (Soifer et 
al. 1987) and the sample "S17" of Saunders et al. 
(1990). 
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Fig. 13. — Plot of the cosmic star formation his- 
tory. The filled triangles and squares represent the 
SFRs derived from the Ha (Gallego et al. 1995, 
Tresse & Maddox 1998, Yan et al. 1999) and UV 
continuum (Lilly et al. 1996, Madau et al. 1996, 
Connolly et al. 1997, Steidel et al. 1999) luminos- 
ity density, respectively. All UV density derived 
values are scaled up by a factor of three to correct 
for extinction (see Steidel et al. 1999). The large 
empty square represents the SFR estimated for the 
faint submm population at z = 1 ~ 3 (Barger et 
al. 1999). The large filled circle represents the 
local star formation density derived from the inte- 
grated FIR luminosity density, which agrees very 
well with the SFR derived from Ha measurement 
by Gallego et al. (1995). Two representative mod- 
els of the star formation evolution by Blain et al. 
(1999; dashed hne) and by Tan et al. (1999; solid 
line) are also shown after a re-normalization to the 
new estimate of the current star formation rate. 
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Fig. 14. — A comparison of the cumulative lumi- 
nosity density for the 1.4 GHz radio luminosity 
functions for the IRAS 2 Jy sample with those of 
the "starbursts & spirals" and "monsters" in the 
UGC sample studied by Condon (1989). While ra- 
dio AGNs with lower luminosity are present, they 
do not contribute significantly to the integrated 
radio luminosity density below Lxaghz ~ 10'^'* W 
Hz-i. 
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Fig. 15. — Histogram of 1.4 GHz radio luminos- 
ity distribution for the UGC galaxies studied by 
Condon (1989). The "starbursts" and "monsters" 
show a relatively clean bimodal separation with a 
dividing line near Li^ghz ~ 10^^ W Hz~^. 
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Table 1 

Infrared Luminous Galaxies with log Lfir > 11.3L0 



Name 


< cz> 


m 


log Li.iGHz 


log ieOfim 




q 


notes 




km/s 


mag 


(W 










IRAS 00057+4021 


13516 


16.8 


22.47 


11.19 


2, 


,872 




IRAS 00091-0738 


35509 


18.4 


23.13 


11.79 


2, 


.807 


jet-like extension to south 


IRAS 00188-0856 


38550 


18.4 


23.74 


11.88 


2, 


.315 


interacting 


IRAS 00234-0312 


20219 


17.0 


23.17 


11.40 


2, 


.409 


VV 837, pair 


IRAS 00262+4251 


29129 


16.0 


23.72 


11.66 


2, 


.079 


LINER 


IRAS 00267+3016 


15119 


14.8 


23.24 


11.28 


2, 


.230 


MRK 551, blue compact 


IRAS 00308-2238 


18381 


17.0 


23.06 


11.12 


2, 


,325 




IRAS 00335-2732 


20771 


17.5 


23.02 


11.54 


2, 


,633 


starburst 


IRAS 00456-2904 


33060 


17.5 


23.71 


11.70 


2, 


,185 


interacting, starburst 


IRAS 00537+1337 


24673 


16.5 


23.37 


11.35 


2, 


,206 




IRAS 01003-2238 


35286 


18.9 


23.57 


11.70 


2, 


,258 


HII/WR, Sy2 


IRAS 01053-1746 


6016 


14.3 


23.29 


11.19 


2, 


,084 


IC1623, pair, HII 


IRAS 01106+7322 


13217 


- 


22.98 


11.16 


2, 


,503 


Spiral 


IRAS 01160-0029 


14274 


15.1 


23.01 


11.11 


2, 


,306 


galaxy-pair 


IRAS 01173+1405 


9362 


14.9 


22.97 


11.25 


2, 


,416 


HII 


IRAS 01249-0848 


14592 


15.5 


23.23 


11.12 


2, 


,142 


MRK 995 


IRAS 01257+5015 


15522 


- 


23.09 


11.02 


2, 


,220 


Spiral 


IRAS 01264+3302 


24145 


- 


23.06 


11.36 


2, 


,489 




IRAS 01364-1042 


14250 


- 


22.84 


11.22 


2, 


,703 


LINER, 


IRAS 01418+1651 


8245 


16.3 


22.77 


11.22 


2, 


,599 


pair, HII+LINER 


IRAS 01484+2220 


9705 


13.7 


23.18 


11.12 


2, 


,185 


NGC 695, pec 


IRAS 01572+0009 


48869 


15.2 


24.14 


12.00 


2, 


,012 


Syl 


IRAS 02153+2636 


15016 


- 


23.04 


11.05 


2, 


,268 




IRAS 02203+3158 


10142 


14.9 


23.06 


11.09 


2, 


,286 


MRK 1034, pair of Syl 


IRAS 02290+2533 


15557 


- 


23.25 


11.06 


2, 


,080 




IRAS 02378+3829 


15052 


15.7 


22.94 


11.11 


2, 


,369 


pair 


IRAS 02483+4302 


15571 


- 


23.02 


11.29 


2, 


,502 


QSO-galaxy pair 


IRAS 02597+4455 


16087 


- 


22.94 


11.18 


2, 


,384 


Spiral 


IRAS 03119+1448 


23006 


17.3 


22.90 


11.33 


2, 


,597 




IRAS 03158+4227 


40288 


17.7 


23.68 


12.11 


2, 


,587 




IRAS 03503-0328 


21730 


- 


23.09 


11.44 


2, 


,520 




IRAS 03521+0028 


45622 


18.0 


23.48 


11.99 


2, 


,720 


LINER 


IRAS 03593-3149 


21462 


- 


22.99 


11.29 


2, 


,497 




IRAS 04154+1755 


16659 


- 


23.61 


11.31 


1, 


,897 


one-sided jet, Sy2 


IRAS 04232+1436 


23972 


- 


23.56 


11.55 


2, 


,173 


LINER, 


IRAS 05083+7936 


16288 


15.8 


23.38 


11.42 


2, 


,280 


patchy, compact 


IRAS 05189-2524 


12760 


15.5 


23.01 


11.59 


2, 


,718 


Sy2, pec 


IRAS 05246+0103 


29037 


17.0 


23.48 


11.59 


2, 


,232 


GPS 


IRAS 06076-2139 


11226 


- 


22.77 


11.16 


2, 


,583 




IRAS 06102-2949 


18297 


- 


23.44 


11.43 


2, 


,214 


triplet in cluster 


IRAS 06206-3646 


32390 


- 


23.87 


11.62 


1 


,919 




IRAS 07045+4705 


19636 


- 


23.27 


11.23 


2 


,14,'; 




IRAS 08030+5243 


25006 


17.0 


23.33 


11.53 


2 


,398 




IRAS 08071+0509 


15650 


- 


23.28 


11.32 


2, 


,244 


radio jet, HII 


IRAS 08323+3003 


17885 


18.0 


22.74 


11.26 


2, 


,693 




IRAS 08340+1550 


23414 




23.05 


11.31 


2, 


,459 




IRAS 08344+5105 


29028 




23.20 


11.49 


2, 


,478 




IRAS 08507+3520 


16748 


15.0 


23.38 


11.03 


1, 


,925 


UGC 4653, triple 


IRAS 08572+3915 


17480 




22.50 


11.61 


3 


,218 


LINER, Sy2 


IRAS 08573-0424 


18425 


17.0 


22.94 


11.19 


2, 


,441 




IRAS 08579+3447 


19645 


16.0 


23.31 


11.29 


2, 


,210 


Spiral 


IRAS 09014+0139 


16081 


17.2 


23.06 


11.14 


2, 


,298 




IRAS 09018+1447 


14847 




23.69 


11.25 


1 


,788 


Mrk 1224, quadruple 


IRAS 09061-1248 


22073 


15.8 


23.30 


11.51 


2, 


,417 




IRAS 09111-1007 


16438 


16.1 


23.39 


11.56 


2, 


,369 


radio double, interacting 


IRAS 09126+4432 


11773 


14.9 


23.04 


11.17 


2, 


,364 


UGC 4881, pair, HII 


IRAS 09192+2124 


23258 




23.06 


11.42 


2, 


,516 




IRAS 09320+6134 


11762 


15.5 


23.70 


11.48 


2, 


,007 


UGC 5101, merger 


IRAS 09432+1910 


16048 




22.96 


11.24 


2, 


,452 




IRAS 09433-1531 


15808 


16.2 


23.17 


11.31 


2, 


,296 




IRAS 09583+4714 


25717 


16.3 


23.72 


11.52 


1 


,959 


Syl+Sy2 


IRAS 10035+4852 


19421 




23.36 


11.51 


2, 


,340 


pair 


IRAS 10039-3338 


10223 


15.0 


22.75 


11.28 


2, 


,667 


IC 2545, pair 


IRAS 10173+0828 


14390 


17.5 


22.69 


11.36 


2, 


,808 


megamaser 
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Table 1 — Continued 





Name 


< cz > 
kin/s 


m 
mag 


log Ll,4GHz 

(W Hz-^) 


log L60M"^ 


q 


notes 


IR,AS 


1 n 1 nn i i ooo 


zzys ( 


17.5 


23.31 


11.50 


2.413 






lUoi i-|~oOU ( 




15.0 


23.58 


1 1 .o4 


2.007 


compact 




1 n'^TG_i_i 1 no 
iUv> / o-|-i iuy 


40843 


17.0 


23.51 


1 1 .82 


Z.444 




IR-AS 


IH/inC 1 A A '~) A 

iU4yo-|-44^4 


27616 


17.0 


23.56 


11.66 


2.331 




IRAS 


iUoUZ — io4o 


16067 




23.01 


11.29 


2.470 


Sy2 


IRAS 


iUobo+244o 


12926 


16.0 


23.31 


11.56 


2.427 


HII, interacting 


TT3 A Q 
IrtAo 


1 1 m n_i_/i 1 n? 


1 noKH 

lUoOU 




zz.yo 


1 1 nsa 
11. uo 


z.o ( y 


Arp 148, ring galaxy 


TT? A Q 


1 1 m 1 Q/i 1^ 


1 /I QK1 
14oOl 




ZO.OZ 


1 1 OO 

1 l.ZZ 


1 Qfl/I 

i.yu4 


pair 


IRAS 


1 1 nCO 1 T71 1 

i iuoy+z / i i 


21080 


17.0 


23.38 


1 1.23 


2.109 




IRAS 


1 1 nn m q o 
i iUyo — Uzoo 


31796 


17.4 


23.75 


11.77 


2. 154 


merger 


IRAS 


i izo 1 -|-oooU 


3101 


12.4 


23. 14 


1 1.32 


2.327 


iNLiO ooyOj merger 


IRAS 


I LA / o — UoU / 


15687 


14.9 


23.02 


11.07 


2.283 






i io iV-\-oooo 


27648 




23.46 


1 1 .46 


2.202 




IRAS 


1 1 1 1 Q Q 1 
1 iOUD+iooi 


o8zUd 


17.1 


23.64 


11.81 


2.355 




IRAS 


iiozy+oUoU 


13070 


16.8 


23.05 


11.09 


2.269 


UtjO D8yt), quartet 


IRAS 


il5o4H-iu4o 


19364 




22.74 


11.22 


2.666 




IRAS 


lloyo — U112 


45177 


17.9 


23.82 


11.95 


2.294 


X-ray source? 


IRAS 


izU / i — U444 


38480 


17.8 


23.42 


1 1.84 


2.551 


G-.rO 

by2 


IRAS 


TOT in 1 r\or\n 

izi iz-|-uouo 


21703 


16.9 


23.38 


11.85 


2.645 




IRAS 


1 o 1 on 1 c o Q o 
izizU+Dooo 


18216 


15.4 


23.21 


11.25 


2.305 




IRAS 


izlYo — oo4i 


17007 




24.24 


11.07 


1.081 


interacting pair 


IRAS 


izioo-|-i io4 


20716 


17.1 


23.36 


1 1 .34 


2.196 


galaxy-pair 


IRAS 


izzoo-|-uziy 


47469 


13. 1 


27,42 


11.95 


— 1.277 


ocz / o, Dlazai 


IRAS 


1 o c /I n 1 c vno 


12518 


14.1 


24.02 


11.98 


2.099 


iVlHls. zol, oyi 


IT) A Q 


loll l-|-4o4o 


1 / zyu 




OQ on 


1 1 oi 
11. Zl 


O OA 1 
Z.Z41 




I'D A Q 


1 Q 1 Oft 1 ftOOO 


9313 


15.1 


23.03 


11.25 


2.384 


JlC^C^ QQOCC -nnl-n XJTT 

ULiLj oooo, pair, rill 


IRAS 


loloo+v>4zo 


6894 


14.8 


23.03 


11.11 


2.302 


11.^ ooo, merger, LillNrljrt 


IRAS 


1 Q no 1 /I 
iozzo — ZDi4 


18468 


16.2 


22.77 


11.15 


2.556 


merger? starburst 


IRAS 


iOvioo— i f UU 


14990 




22.88 


11.18 


2.460 


starburst 


IRAS 


1 Q Q Q c nn /I 


17859 




23.19 


11.16 


2. 194 


U(-jU ooo4, triple 


IRAS 


io4zo+oDUo 


11180 


15.0 


23.59 


11.73 


2.282 


iVmiv ZfO, Oyz, IjliNJjjrt 


TT? A Q 


lo4oo — UUo ( 


23058 




QQ OA 


1 1 QO 
1 1.OZ 


o onn 




TT3 A Q 


1 1000 

lo4ol + l^oz 


36341 


17.0 


26.18 


11.68 


— 0.341 


oy2, merger 


I'D A C 


lo4t)4 — ZtiOO 




18.6 


23. 11 


1 1.34 


2.450 


T TTSTTPD 

LiliN ll/rt 


I'D A Q 
IrtAo 


1/11 riyi 1 Q ten 
141U4 — loOU 


ZzU4U 


16.1 


23.38 


11.37 


2.195 




IRAS 


14147 — zz4o 


23850 




23.15 


11.41 


2.433 




IRAS 


14iOo+z / 4i 


onnno 

zuyoz 


16.0 


23.37 


1 1 .26 


2.118 




IRAS 


1 /I nnn o^ron 

i4zyu— z f zy 


17225 




22.87 


11.09 


2.451 




IRAS 


I/IQ/IO 1/1/1*7 

14o4o — 144 f 


24677 


16.6 


23.67 


11.87 


2.368 


LINER, merger 


IRAS 


14ool — iyo4 


O O 

zD88b 




23.36 


11.48 


2.325 




IRAS 


1 /I OTO Oft C 1 


20277 


16.3 


23.49 


11.68 


2.354 


E? 


I'D A Q 
IrtAo 




31450 


17.2 


23.96 


11.57 


1.776 




TD A O 

IrtAo 


14o4 / -|-^44o 


10166 


14.6 


23.32 


11.08 


2.056 


ULrO yoi8, galaxy pair 


I'D A G 
IrtAo 


10Ulo+z41 * 


20685 


16.5 


23.26 


11.36 


2.299 




IRAS 


1o1Do+4zOO 


12049 


14.9 


23.20 


11.39 


2.353 


iVlrtlv o4o, rlll-|-ljli>lr!jrt 


IRAS 




14172 


16.0 


22.32 


11.17 


2.984 




IRAS 


15233+0533 


16227 




22.86 


11.24 


2.535 




IRAS 


15238-3058 


14318 




23.28 


11.06 


2.095 




IRAS 


15245+1019 


22687 


15.8 


23.27 


11.52 


2.452 




IRAS 


15250+3609 


16602 


16.0 


22.95 


11.58 


2.752 


ring galaxy, LINER 


IRAS 


15320-2601 


20549 




22.98 


11.40 


2.585 




IRAS 


15327+2340 


5534 


14.4 


23.33 


11.77 


2.607 


Arp 220, Sy2, LINER 


IRAS 


15462-0450 


30150 


16.4 


23.41 


11.72 


2.459 


disturbed pair, Syl 


IRAS 


15473-0520 


16994 




22.51 


11.17 


2.801 




IRAS 


16090-0139 


40029 


16.6 


23.87 


12.03 


2.344 


LINER 


IRAS 


16133+2107 


27206 




23.18 


11.48 


2.482 




IRAS 


16161+4015 


23288 




23.17 


11.36 


2.405 




IRAS 


16305+4823 


26329 




23.04 


11.40 


2.523 




IRAS 


16429+1947 


24735 




23.25 


11.43 


2.391 




IRAS 


16474+3430 


33418 


16.5 


23.44 


11.66 


2.403 




IRAS 


16487+5447 


31293 


16.5 


23.61 


11.70 


2.264 




IRAS 


16504+0228 


7298 


14.7 


23.69 


11.33 


1.822 


NGC 6240, LINER, Sy2 


IRAS 


17028+5817 


31779 




23.56 


11.66 


2.318 




IRAS 


17132+5313 


15270 




23.18 


11.42 


2.426 


merging pair 
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Table 1 — Continued 



JN ame 


< cz > 




log ljl,4GHz 


log LQo^rn 




q 


notes 




km/s 


mag 


(W Hz~^) 










IRAS 17208-0014 


12836 


15.1 


23.46 


12.00 


2, 


,703 


HII 


IRAS 17409+1554 


22042 




23.05 


11.25 


2, 


.385 




IRAS 17465-0339 


22538 




23.30 


11.39 


2, 


.245 




IRAS 17487+5637 


19687 


16.9 


22.97 


11.27 


2, 


.546 




IRAS 17501+6825 


15357 


15.2 


23.12 


11.06 


2, 


.194 




IRAS 17517+6422 


26151 




23.21 


11.40 


2, 


.396 




IRAS 17574+0629 


32860 


16.5 


23.53 


11.63 


2, 


,263 


E 


IRAS 18214-3931 


34005 




23.11 


11.70 


2, 


,755 




IRAS 18293-3413 


5449 




23.16 


11.30 


2, 


,344 


HII, pair 


IRAS 18368+3549 


34825 


16.3 


23.74 


11.69 


2, 


,180 




IRAS 18443+7433 


40395 


17.3 


23.60 


11.80 


2, 


,358 




IRAS 18470+3233 


23626 


16.0 


23.19 


11.63 


2, 


,569 




IRAS 18544-3718 


22012 




22.79 


11.40 


2, 


,806 




IRAS 19096+4502 


18980 




23.22 


11.25 


2, 


,262 


pair, HII+HII 


IRAS 19115-2124 


14608 




23.45 


11.38 


2, 


,156 


interacting pair 


IRAS 19297-0406 


25674 


16.0 


23.62 


11.94 


2, 


,488 


merger, HII 


IRAS 20010-2352 


15249 


16.9 


22.89 


11.19 


2, 


,488 


Sy2 


IRAS 20046-0623 


25362 


16.3 


23.32 


11.59 


2, 


,444 


interacting pair 


IRAS 20087-0308 


31600 


16.1 


23.11 


11.92 


3, 


,009 


merger 


IRAS 20210+1121 


16905 


15.8 


23.52 


11.29 


1, 


,873 


interacting, Sy2 


IRAS 20414-1651 


26107 


17.1 


23.55 


11.76 


2, 


,378 




IRAS 20550+1656 


10900 


15.2 


23.05 


11.44 


2, 


,523 


IIZw096, interacting 


IRAS 20567-1109 


15221 




23.26 


11.16 


2, 


,068 




IRAS 21316-1729 


17050 




22.59 


11.23 


2, 


,812 




IRAS 21396+3623 


29034 


16.5 


23.41 


11.51 


2, 


,328 




IRAS 21442+0007 


22187 




23.10 


11.27 


2, 


,413 




IRAS 21504-0628 


23263 


15.6 


23.28 


11.56 


2, 


,407 




IRAS 22491-1808 


23312 


16.2 


22.84 


11.74 


3, 


,025 


merger, HII 


IRAS 23050+0359 


14205 




22.85 


11.13 


2, 


,474 




IRAS 23204+0601 


16779 


16.0 


23.08 


11.33 


2, 


,407 


blue E 


IRAS 23327+2913 


31981 




23.27 


11.58 


2, 


,493 




IRAS 23365+3604 


19338 


16.3 


23.36 


11.71 


2, 


,514 


disturbed, LINER 


IRAS 23410+0228 


27335 


17.8 


23.02 


11.50 


2, 


,608 


disburbed, X-ray source? 
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Table 2 
Radio-excess Objects 





< cz > 




log L/14QHZ 


log ^60^m 


Q 


notes 




km/s 


mag 


(W Hz-1) 








IRAS 03164+4119 


5486 


12.4 


25.13 


10.57 


-0.463 


3C84, NGC1275 


IRAS 03208-3723 


1780 


9.9 


22.24 


9.26 


1.511 


NGC 1316, Fornax A 


IRAS 05497-0728 


2334 


14.0 


22.54 


9.65 


1.301 


NGC 2110, Sy 2 


IRAS 05511+4625 


6009 


13.9 


23.28 


10.32 


1.254 


UGC3374, Sy 1 


IRAS 06097+7103 


4050 


13.8 


23.59 


10.06 


0.618 


UGC3426, Sy 2 


IRAS 06488+2731 


12270 


14.9 


23.71 


10.83 


1.253 


Sy 2 


IRAS 07224+3003 


5730 




23.01 


10.25 


1.434 


UGC3841. FR-I RG 


IRAS 08400+5023 


3226 


12.4 


22.41 


9.57 


1.549 


NGC2639. Sy 2 


IRAS 10459-2453 


3730 


14.8 


22.39 


9.77 


1.616 


NGC3393, Sy 2 


IRAS 11220+3902 


2080 


12.0 


22.02 


9.19 


1.535 


NGC3665. radio jet 


IRAS 12080+3941 


995 


11.5 


22.19 


9.81 


1.388 


NGC4151, Sy 1.5 


IRAS 12173-3541 


17007 




24.24 


11.07 


1.083 


pair 


IRAS 12265+0219 


47469 


13.1 


27.42 


11.95 


-1.277 


3C273, QSO 


IRAS 13184+0914 


9565 


15.1 


23.30 


10.54 


1.520 


NGC5100, pair 


IRAS 13197-1627 


5152 


14.6 


23.19 


10.43 


1.402 


Sy 2 


IRAS 13443+1439 


6435 


15.5 


23.23 


10.35 


1.379 


Mrk 796, Sy 2 


IRAS 13451+1232 


36341 


17.0 


26.18 


11.68 


-0.341 


4C12.50, Sy 2 


IRAS 13536+1836 


14895 


14.8 


24.26 


10.95 


0.844 


Mrk 463, Sy 1 


IRAS 14106-0258 


1753 


13.6 


22.35 


9.69 


1.506 


NGC5506, Sy 2 


IRAS 15060+3434 


13500 




23.72 


10.93 


1.475 


VV 59, pair 


IRAS 21497-0824 


10330 




23.61 


10.79 


1.363 




IRAS 22045+0959 


7800 


15.1 


23.17 


10.51 


1.575 


NGC7212, Sy 2 


IRAS 23254+0830 


8698 


13.6 


23.55 


10.88 


1.530 


NGC7674, Sy 2 



Table 3 
Infrared-excess Objects 



Name 


< cz > 


III 


log Li.4GJfi 


log LeOfim 


q 


notes 




km/s 


mag 


(W Hz"^) 








IRAS 02244+4146 


5606 


15.0 


21.54 


10.35 


3.075 




IRAS 02330-0934 


1509 


11.0 


20.36 


9.12 


3.085 


NGC988 


IRAS 04332+0209 


3580 


15.6 


21.09 


9.96 


3.070 


UGC3097 


IRAS 04489+1029 


8416 




21.87 


10.47 


3.086 




IRAS 08572+3915 


17480 


14.9 


22.50 


11.61 


3.218 


LINER, Sy2 


IRAS 12243-0036 


2179 


14.2 


21.62 


10.57 


3.075 


NGC4418 


IRAS 19420-1450 


-57 


9.4 


18.45 


7.79 


3.184 


NGC6822 


IRAS 20087-0308 


31600 


16.1 


23.11 


11.92 


3.009 


merger 


IRAS 22491-1808 


23312 


16.2 


22.84 


11.74 


3.025 


merger, HII 
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Table 4 

IRAS 60/im Luminosity Distribution and Luminosity Function 



log Leoum 


N 


log PeOfim 








(Mpc~^ mag~^) 





7.875 


2 


-2. 


8.125 


6 


-2. 


8.375 


16 


-2, 


8.625 


34 


~2, 


8.875 


68 


-2, 


9.125 


95 


-2, 


9.375 


128 


-2. 


9.625 


182 


-3. 


9.875 


203 


-3, 


10.125 


277 


-3, 


10.375 


242 


-4, 


10.625 


210 


-4, 


10.875 


138 


-5, 


11.125 


83 


-5, 


11.375 


57 


-6, 


11.625 


38 


-6, 


11.875 


17 


-7, 


12.125 


4 


-8, 



2.44^ 
2.34^ 
2.291 
2.34j 
2.4lj 
2.64+ 



23 
.53 


0, 


.205 


± 


0, 


.320 


15 
.23 


0, 


.307 


± 


0, 


.226 


10 
.12 


0, 


,449 


± 


0, 


.168 


07 
.09 


0, 


.464 


± 


0, 


.117 


05 
.06 


0, 


.491 


± 


0, 


.085 


04 
.05 


0, 


.490 


± 


0, 


.072 


04 
.04 


0, 


.411 


± 


0, 


.057 


03 
.03 


0, 


.447 


± 





.050 


03 
.03 


0, 


.495 


± 





.049 


.03 
.03 





.506 


± 





.043 


03 
.03 


0, 


.507 


± 


0, 


.046 


03 
.03 


0, 


.547 


± 


0, 


.051 


04 
.04 


0, 


.480 


± 


0, 


.059 


05 
.05 


0, 


.467 


± 


0, 


.075 


05 
.06 


0, 


.534 


± 


0, 


.097 


07 
.08 


0, 


.525 


± 


0, 


.118 


.09 
.12 


0, 


.492 


± 


0, 


.170 


.18 
.30 


0, 


.387 


± 


0, 


.311 



Table 5 

schechter parameters for the iras 60 yum luminosity distribution and function 



Functions 


L' 


a 


p'V 




and Fits 






or p* 




n(L)— low L 


2.24 ± 0.41 


-0.805 ± 0.058 


509 ± 160 


6.2 


n(L)— high L 


21.46 ± 1.82 


-0.805 ± 0.058 


143 ± 8 


4.1 


n(L)— single fit 


15.38 ± 0.74 


-1.221 ± 0.019 


309 ± 35 


232 


p(L) — low L 


2.33 ± 0.48 


-0.821 ± 0.072 


(2.23 ± 0.29) X 10"*' 


4.4 


p(L)— high L 


21.69 ± 1.85 


-0.821 ± 0.072 


(2.21 ± 0.13) X lO"** 


4.1 


p(L)— single fit 


15.38 ± 0.73 


-1.220 ± 0.019 


(0.81 ± 0.07) X 10"^ 


234 
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Table 6 

1.4 GHz Radio Luminosity Distribution and Luminosity Function 



log Li,iGHz 


N 


log Pl.iGHz 


V/Vm 


(W Hz-i) 




(Mpc~^ mag"'^) 





19.625 
19.875 
20.125 
20.375 
20.625 
20.875 
21.125 
21.375 
21.625 
21.875 
22.125 
22.375 
22.625 
22.875 
23.125 
23.375 
23.625 
23.875 
24.125 



10 
35 
49 
66 
94 
128 
186 
228 
294 
241 
181 
115 
71 
46 
32 
4 
3 



-2.30"' 



-2 64+°-^° 

-2A3+_°o\ll 
-2.46lS-t 
n+o.os 

-0.10 

-2 49+""^ 

'^■^^-0.08 

-2 54+"-''^ 

"^■'=-0.06 

-2 96+°"'' 

^■^°-0.05 
o 1 1 +0.04 
"^■-^-^-0.04 
I +0.04 
-0.04 



-3.4i: 



r, cn + O.O.S 



-4.06: 
-4.52+ 



- + 0.04 
-0.05 

r, + 0.06 

-0.08 
-5 10+°'"' 

O.OD_Q^2 

-6 00+° " 

D.UU_Q 20 
_K qc+0.22 

o.ao_o 49 
-7 96+°-^° 

- '■0°_1.03 



0, 


.211 


±0, 


.265 





.284 


±0, 


.188 


0, 


.427 


±0, 


.207 


0, 


.444 


±0, 


.113 


0, 


.474 


±0, 


.098 


0, 


.428 


±0, 


.081 


0, 


.489 


± 0, 


.072 


0, 


.406 


±0, 


.056 


0, 


.486 


±0, 


.051 


0, 


.495 


±0, 


.047 


0, 


.512 


± 0, 


.042 


0, 


.505 


±0, 


.046 


0, 


.535 


±0, 


.054 





.465 


±0, 


.064 


0, 


.534 


±0, 


.087 





.494 


±0, 


.104 





.501 


±0, 


.125 





.753 


±0, 


.434 





.561 


±0, 


.433 



Table 7 

Schechter Parameters for the 1.4 GHz Radio Luminosity Function 



Functions 


L' 

(10^^ W Hz-^) 


a 


P* 




p(L) — low L 


2.11 ± 0.28 


-0.633 ± 0.051 


(3.17 ± 0.17) X 10"* 


5.9 


p(L)-High L 


14.43 ± 1.65 


-0.633 ± 0.051 


(0.83 ± 0.08) X 10-*^ 


7.5 


p(Z/)— single fit 


10.76 ± 0.84 


-0.988 ± 0.028 


(2.43 ±0.08) X 10-^ 
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